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13 aBsTmacr 


A threshold detection system based on the Neyman- 
Pearson criterion is derived for sn infrared nutating 
optical system. Detection is optimum ror small signal- 
to-noise ratios and a Gaussian uncertainty in the pointing 
error of the optical system. The signal spectrum and back- 
ground power spectral density for a rectangular space 
filter are computed numerically and used ta specify the 
matched filter for the threshold detection system and 
evaluate its performance. 
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1, INTRODUCTION 


An infrared (IR) detection system determines the 
presence or absence of a target in a particular area of 
space by processing the received IR cnergy from thut area, 
The target must be found amid ever-present background 
radiance and any internal noise generated by the detection 
System, In this thesis, one type of system, called a 
nutating optical system, jis considered, 

Infrared detection is a combination of spatial frequency 
processing and temporal processing. Spatial frequency process- 
ing is a means of achieving a modulated output from the IR 
radiance passing into the optical system so that information 
may be eusily extracted (Ref. 14]. To extract this informa- 
tion, some form of temporal processing is used. Ina 
nutating optical system, the optical axis rotates about an 
axis perpendicular to the image plane in such a way that a 
point image traces out a circle in the image plane, The 
Spatial frequency filter is a piece of IR sensitive material 
in the image plane which produces a voltage output that is 
amplified and fed into a temporal filter for processing. 

A nutating optical system is frequently used in a missile 
seeker because it produces a signal when there is no 
tracking crror and the rotating optics can be made into a 
gyro that provides the inertial reference, 

This work is a synthesis and extension of two recently 


published articles in othe literature, The first, by 





ee SS a 


Harger [Ret. 3], fomiuiares the detection theory of a 
known signal in background ana white noise. Since in IR 
target detection the signal is never known completely, 
Harger's work is extended to che case where the target 
amplitude and position are unknown, However, a known 
target shape is assumed as well as a priori statistics of 
amplitude and position, Since detection is most difficult 
when the amplitude is small, a test (called the threshold 
detector) is derived that is optimum for weak signals. 

To specify the form of the threshold detector and to 
Calculate ity performance, an cigenvalue integral equation 
Must be solved. Ordinarily, this is a difficult undertaking 


but the solution becomes trivial if tho cevariance function 





of the background process is periodic. This is exactly = 
the condition realized for the background process out of 
a nutating detector provided that the phase, which is 
physically unimportant, is averaged out. 

Samuelsson's results [Ref. 11] specify the form of the 
Signal and background coefficients in terms of the system 
parameters and signal and background radiance, These 
derivations were first published in a Swedish interna] 
report (Ref. 10], then re-derived in English by Yoshitani 
[Ref. 6]. The latter's derivations are included in the 
appendices. 

Samuelsson's cquations have been applied to a rectangu- 
lar IR detector assuming a Wiener spectrum for the random 


background and the coefficients heve been calculated, 


a eile a He Rebels cimemarin” = omen oe = 
: ae wae Rages 3 Pdieaineke 
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Using these cocfficients, the optimun filter has been 


Sten 


determined and the probability of detection computcd as 
} ) 


a function of signal-to-noise ratios and background-to- 


noise ratios, 
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‘ IJ]. SPATIAL PROCESSING 


The spatial frequency filtering portion of the detection 
: system is comprised of a lens system and an IR sensitive 


detectei located in thc image plane, the focal plane 


of the optics. To adequately derive information 
from the IR rediation, the Jens system and detector must 
be of a nature that causes the veltage output to be modulated 


in a way that information about a target within the field 


eames Sonematiice mt ci sac. antimssteais e d + 


of view can be processed. In @ practical nutating optical 
system, tne lens system moves and the detector is stationary. 
However, for conceptual purposes, this is equivalent to a 


stationary lens systen and nutscting detector. 





A. TARGET SPECTRUM 


The object to be detected, the target, is located 


or : ‘i 


jin the object plane a distance R from the lens and parallel wi 


to the image plane as shown in Figure 1.1. Tie angle rrom 


the perpendicular axis connecting the planes to the target 





ccordinates measured in the x-direction is called x and is pe 
measured in milliradcians; the angle in the y-direction is 
y. The target is considered to be ciose to the optical 
axis and at a sufficient distance, R, such that small angle 


approximations may be made, 


e (2.1) 
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Consider a target with a radiance distribution 
N(r) (W/em?-sr) which is the radiant powcr into a unit solid 
angle per unit projected area of the sources. r =r(x,y) is 
a vector whose origin is the perpendicular axis connecting 
the coordinates. 
| Assuming no atmospheric attenuation (this can be 
included in N(r) if it exists) and perfect optics, the power 


distribution in the image plane is 


Po(r) + Nr) + Ag Gv/sr) (2.2) 


where Ao is tne effective entrance area of the optics. 

Dividing by Aga one can see that the radiance 
distribution in the imape plane is numerically cquual to the 
radiance distribution in the cbject plane. 


However, because the optics is not perfect, a 


point cbject given by 
Hp 6(r) = dp d(x} Sfy, () , (2.3) 


? 
(where Hp (W/cm™) is the irradiance at the optics received 


from the point source) is imaged as a power distribution 


Hp + £0 (x) (W/cm? Sr) (2.4) 


zs igs the point spread function of the 


where f(r) (sr- 
optics. Irradiance is defined as the radiant flux indicent 
on a surface of unit area [Ref. 5]. 

Therefore, the radiance distribution, N(r), in the 


object plane will be imaged as 


N'(r) = [sors fy (s) d’s (Ww/em?Sr) (2.5) 
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where d?s = dxdy and the integration is over the entire 
image piane. 
The detector is a space filter in the image plane 


which may be a moving retical or simply an aperture across 


a ee ee Se cate weiner ser 


which the power is integrated. 


The powcr in the image plane incident on the detec- 


tor <1s 
: 2 ate 
H~ fx'(r) w(x) ar (K/cm*) (2.6) 


where t(r) = t(x,y) is the transmittance of the detector. 
When the detector moves as a function of time, the 


power incident on the detector becomes # function of time, 


) ar (W/cm?) (2.7) 


ee mem ee 


H(t) = ys (rT) ire 


image plane. 


A nutating detector moves circularly in the image 
plane dDut its orientation remains fixed, i.e., the coordin te 


system (x',y’) shown in Fig. 1.2 does not rotate. The nut.- 


2S eae een ew ew epee Soe 


tion rad.us measured with respect to the optical axis is 


given by p. 


| : where o(t) describes the movement of the detector in the 


i 

i 

| Assuming a nutating optics and no relative motion 
between the target and the optics, the radiance on the 
detector will be periodic and thus can be expanded in a 


Fourier series 






H(t) = EH e QO<t<T (2.8) 
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where 


T ~jnwot 
H(t) ¢ dt 
QO 


= 
1 
fe 


Replacing H(t) by equation (2.7), 
T -jnwt 
Pa i : ° 2 ! . 
Hy T J dt e d r N (r) (r-p(t)) (2.9) 
The motion describing function, p(t), for circular 

nutation is 

p(t) = p(cosut, sin wot) , (2.10) 
where W, = 2n/T (Crad/sec) is the radian frequency of nutation 
and T is the period for one nutation, 


Substituting (2.10} in (2.9), He can be shown (see 


Appendix A) to be 


Hh = i Jricg) t*(k) € JN J, (270 Vk, e ) d°k (2 11) 
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Figure 1.2. Circular Nutation in the Image Plane. 
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where k = (ky ks two-dimensional spacial Fourier 
transform coordinates [Ref. 8], 
N'(k) =» transform of radiance distribution in the 
image plane, 
tT®(K) = complex conjugate of the transform of the 
optical transmittance function in the detector 
coordinate system, 
¢ = tan! ky ky 


p = nutation radius, 


and JC) is the uth order Bessel function of the first 
kind. It has been assumed that the image plane is very 
large compared to the radiance distribution [from a target 
so that the image plane may be considered infinite for 


mathematical purposes. 


B. CORRELATION OF BACKGROUND NOISE 
The background power incident on the detector is fron 
a sample background scene where a scene is a two-dimensional 


random process characterized by the radius vector r = (x,y) 


ee ee 


from the optical axis. 


tae oe 


Let N(r)} be a sample scene radiance distribution on the 


object plane from an ensemble of scenes which has been given a 


suitable probability structure. The background power incident 


on the detector is | ¥ 
B(t) = JN'(r) t(r-p(t)) e?r . (2.12) | 


the covariance function is defined by taking the expected 


value 





EC((B(t)} - E(B(t)}) (B(u) - ECB(u))) = Rp (tu) (2a os) 


where E denotes the expectation over the ensemble of scenes. 
Because the detector is nutating, the randem process, B(t), 
wil] not be stationery. However, if the detector and 
background do net move with respect to each other, then 


Rg (t,>.t2) is doubly periodic: 


R(t, + MT, t. + NT) = Rg(tysto) (2.14) 


where i aul N are integers. 
It can be shown [Ref. 7) that the form of the covari- 
ance of a doubl,’ period process is 
e jmwoty 


Re(ty ste. 2 Bee ge RDG hes (2.15) 
m,n 


Because the nutation is periodic, over the ensemble of 
scenes, there is little significance in the starting time t. 
Therefore, considering t as a random variable that is 
uniform over the nutation period, one can “average it out" 


of tie covariance function as 


a a ae 


T 
R(t) * t po Rete des (2.16) 


oO 


Substituting (2.16) in (2.15), 


Reyes eB e xi e dt. (2.17) 
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Using the relation 


1 T J(u-vjupt 
buv sw fe dt , (25:18) 
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the correlation function becomes 


4 RWgt 
62ue ; C2s0) 


Ry (1) 7 n 


31 8B 


Equation (2.19) shows that B(t) can be considered a 
wide scnse staticnary process and that it is periodic in 


the mean square sense; i.e., 


Ry (1) = Ry (T+) for any 7 (2.20) 


The power spectral density is the temporal Fourier 


transform of Rp (t). Taking che transform one has 


Spe) + en 


3m 8 


Bo § Cu-nug). (2.21) 


The coefficients, 84> can be related to the optical 


system parameters by 


Ye eon ee 
*ky 


8," Ga ieyceres Ws Ck) 1 (209° k, ) d°k (2522) 


where F, (k) is the transform o. the point spread function, 
and W, Ck) is the Wiener spectrum, the transform cf the back- 
ground correlation function. This derivation was first made 
by Samuclsson (Ref. 10] in an internal Swedish report then 
again by Yoshitani (Ref. 6]. Yoshanti'ts derivation has 


been included as Appendix B, 





= Se eee 


Temm rma. 


Equation (2.19) implies that the background, B(t), can 
be represented by a Fouricr series with uncorrelated 
cocfficients [Ref., 7]. Thus, 

jNuyt 


B(t) = = b, c (2,23) 
n 2 


where the convergence is in the mean square sense. 


The coefficients, bye satisfy 


E(B(t)) n= 0 
E(b_) = : (2.24) 
t n ¥ 0 
Morcover, they are uncorrelated 
BF mn 
E(b_b.) © ; : (2.25) 
0 myn 


C. WHITE NOISE COMPONENT 

In an actual system there will be a certain amount of 
internal noise generated in the detector and preamplifier 
which is additive to any tackground noise or signal. In 
the frequencies of interest the primary source of this 
noise is thermal agitation. This type of noise, calied 


Johnson noise, is assumed to be white and Gaussian. 


Let W(t) be Gaussian white noise with power spectral 


density aca: The power spectrum at the output of the 


preamp with bandwidth B is shewn in Figure 1.3. 





Mag. 
(V/Hz) 


B Frequency 


(Hz) 





Figure 1.3. Power Spectrum of White Noise, 


The RMS voltage measured at the output of the preamp is 


na 


Verms = - 2B = NB (V4) . (2.20) 
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TIT. TEMPORAL PROCESSING 


The output voltage from the detector is temporully pro- 
cessed to determine the presence of a target. The type of 
processor is determined by applying statistical detection 
theory to the known behavior of the nutating system and var- 
jous realistic assumptions about the system and the target. : 
An optimum processor in the sense of small signal amplitudes, 
called a threshold detector, is developed in this section. 
Harger [Kkef. 3] derived the detection theory for the 
known signal case in the presence of additive background and 
white noise components. In IR detection, the target ampli- 
tude and position are usually not known completely. The 
threshold detector is an extension of Harger's ark which 


accounts for the unknown amplitude and position. 


The form of the optimal detector under the Neyman-Pearson 
criterion will be derived for multiple observations, [t i: 
assumed that the image plme and object plane remain parallel. id 
Stationary with respect to each other during the observati ns. 

The detection of the target under the Neyman-Pearson 
criterion becomes a problem of hypothesis testing. The 
task of the detector is one of choosing between two 
hypotheses, Ho that only noise is present and Hy that in | ; 
addition to the noise there is a target present. The design 
of the detector is one that permits the correct choice of 


hypothesis Hy (a detection) with maximum probability of 


detection, Qd, in a fixed probability, Qiap of choosing 





ey eT 
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Hy when Hp is truce (a false alarm). The structure of the 


Neyman-icarson criterion requires forming a likelihood ratio, 
Ayo» which is compared to a known threshold, Veo specified 


for a given false alarm probability. 1i the likelihood 


ratio exceeds the threshold, hypothesis Hy is assumed true 
and if not, hypothesis Ho is assumed true, 


The logic process 
may be written 


(52.1) 
aks 


where the threshold, Veo is erived from the expression 


Qra 7 ; p(2Z{H,) dz 


(3.2) 
t 


and Z is a sufficient statistic of the received data. 


A. DETECTION OF A KNOWN SIGNAL 


The data for processing is a set of M functions of time 


{2 (t), n=1l,...,M} where each function is the output of the 


preamp during one nutation (0,T). T is the nutation period. 


The output may be represented as 


Z(t) - S(t) + B(t) +*Ni(t)3 O<t<T, nel,...,M (3.3) 


S(t) represents the output duc to a target. The component, 


B(t), represents an additive background and is a sample 


function of a random field (B). The relative size of the 


target in the ficld of view is assumed sufficiently small 


that the background is additive while B(t) is assumed to be 


21 





No (t) represents the additive noise component due to the 
internal noise of the detector-preamp combination. It is a 
sample function of a white Gaussian random process (vw) of 
zero mean and power spectral density x/2(ve/11,). 


: : . aT 
the sample function during cach of the M nututions observed. 


written as 


Hy: (2, (t ) = SyCt ) + Bt) + N(t ) 5 meds... M) 
(3.4) 
Ho? SE Ct) SoBe le RCE one dyy eso 


Denote the likelihood ratio as Ay ({2,)) where the like- 


The two hypotheses for the likelihood ratio may be 
| 
lihood ratio is defined as 

1 

( 


Ay, ({2,)) Ses MEN (3.5) 
10 n p({Z} Hy) : 
The likelihood ratio is complicated by the fact =n and 
z,, are not independent because of a common noise component 


B(t). To simplify the derivation Harger introduced an 


auxillary hypothesis 


H: {2 = Nn ; ne=l,...,M}, (3.6) 


Then the likelihood ratio, Ayo» can be computed using 


the chain rule for likelihood ratios: 


PUCZ FjH,) PC{Z PH.) sp C{Z,} 1H) 
Mo" POULIN) ” PUCZFTH Spe Fp) 


Ay2/fo2: 
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If the sample functions for the background, B, were non- 
random, Z,(t) and 25 (t), i # j, would be mutually indepen- 
dent because of white noisc. Accordingly, if B were fixed, 
one can write the conditional likelihood for the nth 
observation At ({2,) |B). Moreover, the likelihood ratio 
for M observations would be the product of the likelihood 
ratios for cach single observation, 


M 
igs n 
Reg Mayr lag Gls 


The conditioning is then removed by averaging over the 


ensemble of backgrounds 
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Since the noise is white Gaussian, the conditional 


likelihood ratio is that for detecting a known signal 


S78, in white Gaussian noise and has the well known forn 
({Ref. 15], page 253) : 
i 
n : 2 M ot 
= e i ° 
Ay. (2,31 B) CXp | fi o : 2 (t) (S, (t) * Bt) } at 
LT sos (t) + B(t))? at! 
N n=] 0 : { } i 
which may be rewritten as 
Z go, <2 | 
Ay, {2,3 1B) Feexpy py & 5 Z(t) S(t) dt 
n=i-0 ; 
7} M i 7 
Pa nf Si(t) dtp exp a (b) (3.14) 
“n=l 0 
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where 


Da oot NM .J.2 
a(b) = % Ra f @ ty +s. @)y Bide = | f B°(t)dt. (575) 


; ; seal 
To rind Ayg (2d, one must evaluate E,[A 122,318) ] 


where 


M T 
n e 2 ; 
Egta y2((2,}1F)] = exp {% - s z(t) S(t) de 
(3.16) : 
Lt ss 2¢s)ae) {exp a(b)} : 
7 E sf S$ “(t)dt} » E, {exp a . 3 
2 n=1 0 m Bp Gs 


Assum'ng the background is Gaussian the expectation can be 





evaluated by r.presenting the background in a Karhunen- 


Loeve (K-L) expansion 


ls zc 


B(t) = Fb, ¢(t) , O< t<T (S27) 


L 
k 


where ¢,(t) is the eigenf nection coricspendirg to eigenvalue 
k g Pp & 


AY of the homogeneous integra! equation 


T 
Chl: Bs xf Ry (t,u) o(ujdu , Ost <T, (3.18) 
0 


ALP, \ 


ew re cernmnae- + 


Ry(t,u) is the covariance function of background 
Ry (t,u) POTBCE} < My(t)][B(u) - mp Cu) ]} S219) 


where m..(t) is the mean of the background. 
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Then Ep(expa(t)) can be calculated in o straightfcrward 


Manner. After the lengthy calculation shown in Appendix C 





: : M 
in ER {exp a(b)} = - ; Ein (1 + Ny Ay) 
k 0 
2 T M T M 
+ f dt £ (2 .(t) - S.(t)) f du © (2 _(u)-S_ (u))q(t,u) 
No ner 2 0 ma  ™ a 
cot ~ M . . 
ie ; dt {m, (t) my t (Z(t) -S)(t))} h(t) (3,20) : 
n=l] : 
where q(t,u) satisfies the integral cquation 


T 
aH q(t,u) + s Rett aGyuldy Rey. Oe weet: (G28 


and h(t) satisties 


jh nana ne ne ee SIs as al Relea a 


T | 
N 
h(t) + SR (t,t)h(t)dt = m(t) O¢<t< T. (3.22) | 
ou 0 B B - "= | 
q(t,u) is the impulse response of the filter that gives | i 
a i . 
the minimum mean square error estimate PB of the backgroun 
when the mean is zero [Ref. 15]. The additive white noise 
is reduced in amplitude by 1/M. Using the k-L representa- 
tion, one can also write 
oo N d 
og BN . 
q(t,u) = 2 oft) of (4) O <u, t <T (3.23) i 
! oS re 4 
“ot 
a 
Wsing equations (2.9), (3.10), (3.14) and (3.29) and 4 
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observing that 


Ago (2,}) = 4,,UZ, HIS, = 0 (3.24) 


it can be shown (see Appendix C) that the likelihood function 


for hypotheses Hy and Ho is 


M T 
Ayg (20) eexp{r f Z (tia, (t)dt + Kt} C3..253) 
n=l 0 
where 
2 » ~ ©» 
q,, t) “iy (S_(t) ~ ; a S, (u)q(t,u) du} (3. 25a) 
g= 
and 
ie : M oT ah , ITM. M ; 
Lr og (t)du + ff = ES (tjyS,.Ca)q(t ,u)dtdu 
n=lo 2 NM oon=1 221 2 
2 M T 
-= = f S (t)jph(t)dt (3.25b) 
™ n=1 0 o 


The logarithm is a monotonic function. Thus it offers 


no loss of information but greatly simplifies the calcula- 


‘tions. In the remaining work 

A" yo ({2 3) = ln Ato etd 1) i (3.26) 
Rewriting equation (3.25), 

At, )({2,}) = Lb ({2,3) + Kt (3.27) 


where 
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K' does not depend on the input data so it may be included 


in the threshold. The hypothese:, test now becomes 





on comparing q,, 8 and (3.21), it is seen that the above 


n 


model approaches the Gaussian model if 


Hy 
L( 2, ) 2 Threshold . (3,28) 
Ho 
It is also possible, as Harger has shown, to derive 
tests similar to (3.28) without making the Gaussian back- 
ground assumption, In particular, if the background fluctu- 
ates many times over (0, T) such that a(b) o. (3.15) is : 
approximately Gaussian, then the test is i + 
— 
Hi Ds 
M oT s cS 
| £ oOff 2 (tie, (t)dt Z Threshold (3.29) PS 
n=1 0 H eee 
| 0 s ; 
where 
: | ¥ 
T M . 
1 2M . Le. Xe 
g(t) = y (S(t) - S/S Ry(t,u) [GF 2 $,(u) ida ? 
0 ge] Bs. 
Since by Mercer's theorem ¥ 
Rr i a 
Rz (tu) : AG, (tT) OFC), Os ty, ucT (3.30) | = 
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Tne above condition, which implies a weak background, is 
sufficient but by no means necessary. The two models 
approximate each other whenever the second term of qn and 

8, are nearly identical. The point is that while the 
Gaussian background assumption may not fit the reality, 

a processor based on the assumption is probably not far 

from being the optimum. Also, the Gaussian assumptio’ allows 


one to compute the processor performance. 


B. DETECTION OF SIGNAL WITH UNKNOWN PARAMETERS 

The likelihood ratio derived in section A is applicable 
only if the signal is knovn completely. Such is seldom the 
case in IR target detection and in general the signal is a 
function of some unknown parameters. Typically, the ampli- 
tude, a,of the target and its position e. in the image plane 
are the unknown parameters, and the form of the signal can 
be rewritten to include these parameters: 


i = 
Si, = Sy (tia, Te) 


Furthermore, by assuming that only point targets are of 
interest and including the assumption of fixed image and 
object planes over the number of nutations of interest, 
the shape of the waveform will be the same during each 
hutation. However, the amplitude is allowed to vary between 


successive nutations. 
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The signal wuveform may now be written as 


S(t 38,7) ea. £(t,t) (3.31) 


where a, is the signal amplitude and Str) is the signal 
shape for a target at ce 
Rewriting hypothesis Hy under these assumptions 


gives 


H)': (2,(t) © a-f(0,1) + B(t) + Ni(t); n=l,....M 


0<t<T) (3. 32) 


The likelihood ratio Ayo becomes the conditional 
likelihood ratio Ayo (iz dlast) obtained by replacing 
S, {t) by eof (tyr) everywhere in (3.25). 

Assuming a priori statistical knowledge of a and Ts 
and their independence, the likelihood ratio may be 
averaged over the respective density functions to produce 


an unconditional likelihood ratio 
Ayo ({2,)) = [Jo C2 ,2la.to)pladp(t,)dadr,. (3533) 


* . < > 
Instead of averaging as above, one could estimate a and T,> 


by the maximum likelihood procedure and substitute the 





estimated values a and t in the conditional likelihood 
ratio and try to maximize the generalized likelihood ratio 


Fas However, this procedure is not followed 


here because the maximum likelihood estimates are difficult 


Ayo (2 31a, 


to obtain. 


Since the main interest is in detection of the weak 


target, the small signal case is of primary interest. 


Expanding AypC(t dase) in powers of signal amplitudes 


M 
> 3 > 
Myo tty eek g) od Een ta Mpg AA Seto) lang 
* (3.34) 
. , M M 32 s 
v 
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Integrating (3,34) as in (3.33) term by term and 


vo 


retaining only the first non-zero term that depends on the 
data provides a statistic which is optimum for weak signals 


and is called a threshold detector [Ret. 4]. 


aaa 


The evaluation is made easier by the fact 


3 > 3 * 
A,,({Z }la,r,)| = dh hehe. age 
oa) 10 ad day 10 noes? 0° and 


(3°33) 
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The first non-zero term becomes 


of 
a _ + a: + 
Yay In Ay, (Zi) last ee NY : Z(t) f(t, rp)dt 


; pn od T : 
- ; 2 (t)dt ; q(t,u)£(u,r,)du 


T ss 6 
- f £(t,ro)h(t)dt (3.50) 
0 


; > 
Performing the integration over To» 


M . T : T : 
Eo ad -2.(t) [£(t) = 2 ae, a). tha) dul 
o ® 0 


= 2 
Ayo ({2,)) Pilty n 


n=l 


. T ss 
fF ECC) Wet): dt (a232) 
0 


where the bar denotes the averaged quantities. The first 
and last terms are independent of the data and may be 
included in the threshold value. Morecver, if the ampli- 
tudes are identicaliy distributed, as is generally the case, 


the threshold detector performs the test 


Hy 
M oT ; ; 
Woe f. f Z2.(t) alt) dt 2 W (3.38) 
n=l o % . 
H 
0 
where We is the threshold and 
: oe 7 : 
GUE) Se LEC atte) 2a aay) 2 (5:38) 
0 











yr Ee are 


ee eee cam Pe oa ees 


Thus the threshold detector can be mechanized either as a 
correlator or matched filter followed by an accumulator. 
The filter function q(t) may be represented in a series 
oxpansion in terms of the eigenfunctions of Rp (t,u) and the 
cocfficients of the signal expansion. Expanding f(t) in 


terms of the eigenfunctions: 
f(t) = z f,o,(t) . (3.40) 


Substituting (3.40) in (3.39) and using (3.23) 


ms «f 


q(t) = aH 5 oae o, Ct) . (3.41) 
m * “k 


C. EVALUATION OF PERFORMANCE 

The Gaussian assumption was made in deriving the 
threshold statistic; therefore, the threshold statistic 
itself is Gaussian. This means that the probability of 


detection, Qd, and probability of false alarm, are 


Qrar 
uniquely determined by second-order statistics, e.g., 


eee 


the means and variances of W under hypotheses Ho and Hy: 


Recall that the threshold statistic is 


M T 
We Ef f 2 (t) q(t) dt. (3.42) 
n=] 0 


The means and variances of W are: 





mM oT 
E(W]Ho) = £ s q(t) E(B(t) + Ni (t)) dt 
n=1 0 
MOT, 
=» <I fs Q(t) m(t) dt 
n=l] 0 
MoT. : 
R(WIH)) = © fs q(t) Ela, £(t,t9) + Blt) + Ni (t)] dt 
n=l 0 
Me ; 
= E(W[H)) + £ f £(t) q(t) dt 
n=l 0 


VAR (WI ly) = VAR (WIN, ) = VAR 


g(x) Ra(t,u) dudt 


oO‘ 


(3.43) 


Letting py C° 1H; ) be the Gaussian density function of 


W under Hy which is specified by the equations above, then 


Qra = Py (x1HG) dx (3.44) 
t 


and 


Qg * / PY (x1H,) dx (3.45) 
me 


where W, is the threshold voltage. 


Equations (3.43) and (3.44) can be solved to yield 





Qua * O09) (3.46) 
and 
Qa te. (x-4) (3.47) 
where. 
es We E(WiHG] | 
YVAR 


d is the equivalent signal-to-noise ratio 


E(WIH,] 
/VAR 


M te a2 
E- a. £ £0t) g(t) ‘dt 


Se 
, Q(t) f q(u) Ra(t,u) dudt 
0 





ef a2ct) de + M? 
0 


and $. (8) is the complementary error function of the kind 
2/2 


os) se Le av (3.49) 
vin g 





By setting a desired probability of false alarm, x 
| may be determined from equation (3.46). The probability 

| of detection is then determined by solving for d and using 
equation (35.47). 


By substituting (3.41) and (3.42), and using (3.23), 





the numerator of d is 


34 
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: f £(t) att) rN ite 
roa tha(t)dt = £ r (3.50a) 
n=] 2 Q ne Hy Ne ay 


while the denominator of d (see Appendix D) is 


n |” Ea 4 
YVAR = oN t 4 aaa (3.80b) 
k oy * Ay 


Thereforc the effective signal-to-noise ratio is 


ae i; 1,1? see fe i 
oe 1 N ee a 


Knowledge cf d along with (3.46) and (3.47) completely 
specifies the performance cf the threshold detector. 
The equivalent signal-to-noise ratio for a target whose 


amplitude and position are known is easily shown to be 


Saal ig whe: «eR Ef, (2.32 : ic (34 2) 
Mine Oa” jae OR Oe ey 
WN “k 
By letting the background be zero it will now be shown ‘= 


that d' reduces to the equivalent signal-to-noise ratio for 
a known signal in white Gaussian noise which is well known 
(Ref. 4). 

The energy dissipated during one observation in a resis- 
tor of 1 ohm if the sigral at(tor) is the voltage across 


that resistor is 


35 
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ny k*o : 


The voltage resulting from M observations then is 


VED ee Peale tie. (E212 : (3.54) 
wv _ a T . * 3.54 
Mey 8 Mary BI k*‘‘o | 
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Substituting (3.54) in (3.52) with ee ® 0, 
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d' = vE —M = RE mw , 
wv VN 
which snows the familiar result that detectability depends 
only upon the signal energy and the signal-to-noise ratio _ 


increases as the square root of the number of observations. 
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TV. THRESHOLD DETECTOR FOR NUTATING SYSTEM 


In an infrared detection system the output voltage to 
be processed is obtained from a detector-preamplifier com- 
bination. The irradiance on the detector in W/en? is con- 
verted to volts with a linear scale factor. The voliage 


output can be written as 


v(t,T,) « KH (t,t 9) (4.1) 


where H(t 7.) is the irradiance on the detector from a target 
at position te and K is @ scale factor in (v/w/em?). Also 


the covariance function may he written as 
. 
Ryg(t) © ELKB(t+1) KR(t)} = KORg(r), (4.2) 


where R(t) is the irradiance on the detector due to a nartic- 
ular background scene. The output voltage due to this scene 
is 


Vz, * KB(t) (4.3) 


The problem now is to relate the detector-preamplifier 
output to the equations derived for the threshold detector. 

Jt was shown in Part III that an optimum processor for 
weak signals was a threshold detector that performed a 
certain linear operation on the data. The linear filter was 
Specified in terms of the cigenvalues and eigenfunctions of 


the integral equation (3.18%), where the kernal Ry (t ,u) is 





the covariance function for the backprounc. 1 & 
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In Pert JI, the covariance function was noted 


periodic. Rewriting (4.2) in terms of the output voltage 


jkw, (t-u) 


‘ rae 
Ryp(t.u) = kK . 8, e (4.4) 


Suvstituting (4.4) into integral eouation (3.18), the solu- 


tion for tre eigenvalues and eigenfunctions is trivial with 


jkwot 





o,(t) = § 


and 


me: 
My Set (4.5) 


where T is the nutation period. 
The signal coefficients were also noted to be periodic 


and could be expanced in a Fouricr series 


jkw,t 


v(t) = K 7 Hi (mje (4. ) 


Necognizing the similarity hetween equations (4.6) and (3. 1), 
the signal coefficients for a known signal which is expanded 
using the basis set above hecome 


f(r \ = KH) (rg) ¥T (4.7) 


° 
For the threshold detector, the signal is averaged with 
respect to an a priori density function for Foe Averaging 


can be done with the ccefficients to give 
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‘ fy s f P(ro)f, (rg)ar, 

= s r Tt )dt 

; KYT f Hye (ry )p (ry) ar, (4.8) 

: KvT Fy 

Using equation (2.10), the coefficients Hy can be written 

{ as 

} 

~ 1 an fat act at -jno , | [2av? ja? 

| R= 3 fot pG,) NUKE) ve (ke Jn [27 VERY [A k (4.9) 


The transform of the radiance distribution from a target, 


—-=— 


N'(k,¥) is a function of the target radiance, R(k), the 


bere + ; 
5 position of the target, To? and the point spread function of 


the optics. FO CK). Rewriting (4.9), 
é . 2 cae 4 
f= 5" fatanag cea 3, [200 keg eas 


( 
| ; > > -j20k,T, 
! fer, p(t je (4.10) 


Tt TT, Le Te 


If the density function, P(r.) is assumed Gaussian with a 


staiudard deviation Po? then the averaged coefficients hecome 


' 
b Saye (ke+k2 
' s : -<7p (k +k \ i ae | 5 
Bie, je en : CORY a -jn¢ aha Atel | ae 
Ht J Ree t(k)e Jn [270 kyoky ek. G12) 
The Gaussian assumption on Ts is reasonahie hecause rs 
! often represents a pointing error from some designated tar- 
get position. Then Po describes the accuracy of the pointing 
| mechanism. Substituting (4.8) into (3.42), the filter func- 
m4 tion for the threshold detector becomes 
bt 
et jku_t 
i ; eo ae 1 “e o 
, k 1 + + KTR VT 
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Rewriting the equivalent signal-to-noise ratio in terms 


of equations (4.8) and (4.12) one has 





M r 4 
Money om ee ie 1 
d = ———__~—_——__|[£ Ht | a (4.13) 
YKJ2 k jw EEK 
“W/2 
It will be shown in Part V that B, may be represented a3 

2,2 

Bs op 2° Qy 


where op is the amplitude of the background noise coriela- 


- 


tion function, 9? is the instantaneous field-of-view of the 
detector and Q. is an integral function of the optical and 


detector paramet*rs and correlation length of the backyround 


radiance. 


Define the fellcwing terms: 
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Signal-tco-noise ratio (SNK) = ee (4.1424) 
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OpokyT 
Background-to-noise ratio (BNR) = ——— (4.14b) 
YN/2 


Using (4.14a) and (4.14b), d may be rewritten as 


eA ame Rr 


c 5 
= 2 l ] 
Lk 1 + MQ, RNR™ 
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Noise cquivalent irradiance (NEI) may be defined at the 


detector-preamplifier output to be 





K*NEI = JF - 2B. = WSR (4.16) 


where B, is the detector-preamplifier bandwidth. 


The SNR and RNR may now be written in terms of the NFI as 


<a>YNTB 


: = ‘R! 
NET SNR! YSTB 


YM SNR = 








where NTB is number-time-bandwidth product and 


oO 
YM BNR = WET YTB = RYR' YNTE. 


Expressing d in terms of NEI one has 
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© 5 
d © SNR’ VNTR xr [H,|? eee ee (4.17) 
k 1 + Q NTB RNR'* 
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VV. SOLUTION FOR RECTANGULAR DFTFCTOR 


The equations in Part II for the detector-preamplifier 
output voltage and the correlation function of the hackground 


must be solved to obtain the filter function and specify the 


detection system performance. 


The two general equations obtained in Part II are 


juk-t 


Hy (Fy) * Mfrare aE eGR) Oo 
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where Hoge) is the nth Fourier series coefficient for the 


radiation on the spatial filter and 
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where By is the nth coefficient of the hackground correlation 


5 


‘ 


function at the spatial filter. 
A rectangular shaped detector was chosen for the calcu- 
lations because it exhibited the simplest form of the 


equations. For circular nutation the filter is represented 


graphically in Figure 5.1. 
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Figure 5.1. Nutation with Rectangular Detector. 


If the detector function is given by 


x +w/2 
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1 * -w/2 < x! 
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0 otherwise 


then its transform is 
~“Jam(kyx,tky yo) 


t(k_,k,) = S——,—_— 
ee mk k 





sin (nk, w)sin (tk h) 


Y 
so that 
. 2 bose 
2 sin (Tk w) sin” (rk h) 
|< (k)| 2 a ee (6.1) 
T ky qT ky 


The point spread function of the optics is assumed to 


be Gaussian: 





Zoe 


f(r) > —-y . 
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: Its transform is 


-2n?o? (k24k2) 
F_(k) =e yo, (5.2) 


The correlation of the background radiance is ass med to 


be of the form 





2 & “ 
a Op (xy) = cGpe alx|, B\y | (5.3a) 


1 


where a ~ and qo are the correlation lengths in the x- and 


y- directions, respectively. The Wiener spectrum is the 


Fourier transform of ¢p: A 
2a 28 

Wp (k) = oF SS a (5.3hb) . 
~ a" + (21k) B- + 2mk i)” ] 
Substituting (5.1), ($8.2) and ($.3b) into (2.10), the 
radiance from a point target becomes # 
‘5 (k " ) . 22 k-) : 

a -jen rotk or -27 "oO" (ko+k? ie 

Hu (ray os i” fe ce, hee” 56 aed 

n‘'o 5 
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-jn tan” *(ky/k) —; 
e a ano fae+k: dk dk. (5.4) 


Letting 
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equation (5.4) becomes 
ii ip 
n -(uft+ul) 
“ . x Yy 

H T,) iy fe 

T 
0D 
Io 


o littcome Yuet sy u.. | 


; -1 
-jn tan “(u,/u_) Js 
e YS 0, 22 furour duydu, (8.8) 


If z is random, the radiance from a target may be 
averaged with respect to the density function p(t) as shown 
in Part IV. In particular, if T. is Gaussian distributed, 


with standard deviation Py the averaged coefficients become 


22 oy 2oae P 3 
ies nfo VE CRY) 52m Oy tot ky) 
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Ba nee Soe a -] ; 
Sink, sinzk hh jn tan (ky 7k) 
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where 
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' 
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equation (5.6) becomes 
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And finally substituting equations (5.1), ($.2) and 


($.3b) into equation (2.22), the correlation coefficients 





become bs 
~An2g2(k2ek*) sin? (rkow) sin’ (nk h) 
By, - fe ay ee 
(nk,) (ny) 
anp /ko+k* |dk dk 
a aot e/ 
°s + (2tk)' BY + (tk) ‘| nae 2 
fe oN 
(ve) 
Letting 


uy = 2npk, 
uy = énpky 


equation (5.8) becomes 


7 mow W ' 2 ) h uw 2 
op? abo? “uy Wy sin vay x oe sin oes uy 
Pn 1 : a ul! — vu" 
y2o * YtZo 


1 1 2 p =o° Oy tl 
See Jn G uy uy jane (5.9) 


(ao) ou, (fo J ruy 





The integrals in (5.5), (5.7) and (5.9) were difficult 
to evaluate. Although the components of the integrand are 
Telatively simple functions, the argument of the Bessel 
functicn prohibited separating the integral into two one- 
dimensional integrals. No closed form solutions were found 
So the integrals were evaluated numerically. The numerical 
method used is based on work by Pierce [9] who applied 
Gaussian quadrature formulas to two dimensional integration 


by inteyrating over a planar annulus in the (x,y) planc. 





Gaussian quadrature formulas are a means of evaluating 
the integral by summing weighted values of the integrand at 
specific points. Fer the one dimensional case, 
| b b 
! J g(xjdx = f w(x) f(xjdx (5.229) \ 
a a 
where g(x) is the integrand to be integrated and w(x} is a 
weight function for which the specific integration was 
derived. For well behaved functions the integral may be 
evaluated as 
b m Bes 
£f g(x)dx =  . A,f(x.) + error (3.131) ¥ 
: i i ; 
a i#l ae 
where the x;'S are the m zeros of the mth polynomial 
m 1 ¥ . 
P(x) = TM (x-x,) | 
™m i=] pS ; ee ' 


of the set of polynomials mutually orthogonal on the inter- 
val {a,b] with respect to the weight function w(x). The 


weights are given hy 








b 
A, = J wx) 1, (x)dx $13) 
a 


where 1; (x) ® PrOxd/b(x-x, PL (x3 )] is the Lagrange inter- 
polation cocfficient. 

For a weight function w(x) and interval (a,b) the poly- 
nomial P(x) and its zeros x; and weight factors A; need be 
computed only once. For a number of weight functions and 
intervals the set of orthogonai polynomials is known, Stroud 
and Secrest [13] give x,'s and Aj's for a variety of weight 
functions and internals. The degree of the formula deter- 
mines the number of x,'s and A,'s, A 2M-1 degree formula 
will have M points and is exact for polynomial integrands of 


ere. 


Pierce applied Gaussian quadrature integration te two 


| 

dimensions where the solution is of the form | 
Sfg(x,y)dxdy = r x Dep 8 Osi; 
ij | 


The integration is over an annulus in the x,y plane with 
inner radius R and outer radius 1. Rewriting the integral 


in polar coordinates 


1 20 
f rG(r,@)d8dr = rt Y,.Gi(r,,6;) (5.14) 
° ij 


; ij 7a 


Pierce showed the summation could be rewritten as id 3 


4(m+l) m+] 
or ae Ai B G(r; cos; 4T; sind;) : (5.15) 
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where A; ana B, ure the weights for the radial and angular 
directions respectively ard 4M+3 is the degree of accuracy. 
For an arbitrary annulus the formula may be obtained by 


rewriting equation (5.14) 


an 72 
Tet a f(r,6) rdrd9g (S.16) 
0 r 
1] 
and letting 
Yr 
¢o.2 — 
¥ 
2 2n 1 
Is TS f f pf (rsp ,8) dodé. (S617) 
“ oO 


Yr 
l/r, 


The approximate solution can be derived as 


4(m+1) m+ | —- 
I = ee i A;B,f | ryto,(r3-ry], a (5.18) 
where 
A, = 2n/4(mel) 
and 


s f . = 
B i nad TI )/2, We weight 


M+l is the order of the orthogonal polynomial, in this case 
the Legendre polynomial on the interval (0,1) ands; 's are 
the zeros of this polynomial. This type of formula is known 
as a spherical product Gauss formula [12]. 

The annuli were picked at radii corresponding to the 


zeros of Jy) because Jo) is the fastest damping component 
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of the integrands. One additional annulus inside the first 
zero was found necessary to account for the 2 term when 
ao is small] and x (and y) approach zero. The distances 
between annuli outside the first five to eight zeros were 
found not to be critical and thus spaced arbitrarily. 

The Gaussian Quadrature formula used was a Gauss-Legendre 


24-point, formula as listed in [13] thus giving a degree of 


accuracy of 99. The integration is limited to only the first 





quadrant as mentioned below. Thus for a 24-point formula, 
(M+1)? or 576 points per annulus were used to evaluate the 
integral. 

The expression for BY is easily seen to be an even func- 


tion both radially and about both plancr axes, therefore, 


| 
| 


only integration over the first quadrant was necessary. 

' With a little more difficulty, tr) and bi can be shown 
(see Appendix F) to exhibit the same property provided one 
expression is used for even n and another for odd n. By 
calculating only over the first quadrant computer time was ¥ 
considerably shortened. 

In order to solve for coefficients of arbitrary order, 
the Bessel functions of that order must first be obtained. 
The available library subroutines were found inadequate for 
this use for two reasons. First, these routines have an 5 
upper limit of 100 on the order of the Bessel function that 
can be computed and second, the routine had to be called for 


each individual order. 
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The recursive relation 


to 
3 


JO) * BIO) + Iya) (5 19) 


Ne 


provides a basis for generating an array of Ressel functions 
of various orders, However for accuracy Jig (¥) and Jno) 
must be known. 

The most accurate method to numericaily calculate Resse! 
functions of various orders and argiments was found to he a 


uniform asymptotic expansion involving Airy functions [1]: 








aay fA Fa) = a0) 
ose fe nT yey nek 
At(n?/3,) © by (A) 
ope! Lo (5.420) 
no/3 k*0 nn 
where A; and A} are two types of Airy functions and 
a 
: a3/? in| 1 + viet |/s| - Ni-2° Seri 
a, (4) se | 
by 0) x -5/48 ee bls 2. < i 
VE ee Clee jee BChea sy? 
ay (a) << a, (a) 
k > 1 (5.21) 


by (4) << bo) 
The Airy functions are calculated using the equations 


oo 8h bee 
A(x) © 5 x ef £(-8) 








where 


x3/2 


a 
tt 
ba] ho 


(5.22) 


and 


x= n/3 4 


The functions f(-6) and g(-6) are tabulated and linear 
interpolation is used to determine their values. A is guar- 
anteed to always be positive if n is selected to be larger 
than the argument. 

By picking an order much higher than the argument, v, 
and Joep using the asymptotic expansion, the recursive 
relationship (5.19) may be used to generate an array of 
successive orders of the Bessel function down to and includ- 
ing Jo: However, errer buiids up rapidly using this method. 


A second relation 

1 = Joo) + 2Jo(y) + 25, (7) Pog kid (5.23) 
may be used to generate a normalizing factor 

C827 (3g). 20,09) 3 2g) ves | (5.24) 


Multiplying each Bessel function value generated in the 
recursive relation by C resulted in very accurate values 
being obtained. The computer program used ta solve the 
equations is listed as Appendix H. 

To evaluate the numerical results a check was used which 
summed the coefficients and compr.cd this summation to an 
analytic expression for a sum that could easily be solved 


numerically. 


For a particuler phase angle of the nutation cycle, wt, 


the radiance function for a target is 


e 7 jku ot 


H(wot) . : Hy (rae (S.25) 





For e point target and a rectangular, circularly nutat- 


ing detector 





-an?y*k? sintwk, jak, (ecoswot (tek) ) 
H(wot) = fe ak, e dk, 
-2n207k? sinthk. j2rk (psinu,t- (ry-y,)) 
‘fe y e Y . y 0°" ak 
T 
Y 
(5.26) 
which has the solution (see Appendia F) | 
= 1 : | 
H(w ot) a [Fy (rego . Fy Cry xr] [E, (ry s¥ Qh) : 
: F2(ty 9 it)] | 
where 
r-Z-pcosuyt 4 i 
E, (7,2,a) = Ert|——_— = | 
¥20 2¥20 4 
r-z-psinw t s 
E,(1,z,4) = Frf a + 2 ; 
4260 2v% J 
| and Erf is the error function ; 
2 : “is 
Erf(x) = —-f e du . (S,27) 
Yn 0 


By evaluating the analytic expression at the phase in- 


stant (wt) that the detector crosses the target and comparing 


; ; ( 





a 





it against equation (5.27) for ar appreciably high n, a 
reasonably check was made on the accuracy of the individual 
Signal coetficients. 

For the coefficients averazed over te the check was 


made at wot = 0 where 
H(o) = 5H, (5.28) 
kee = A 


and the anaiytic expression is 
H(o) = F [Fy Cry okgoW) ~ By (ry exy ow) [Fs (zy 7 gh) 
: Es (1, s¥qsh)] 
where 
E,(r,z,0) = F2(T 1209) |y exo (5.29) 


The rms background voltage was used to check the coef- 
ficients in the correlation function. The mean square volt- 


age is 


ef 
B 


wD 


= E(B2(t)} = Ra(2) = (5.30) 


we 8 
Ye 


Appendix B shows that the correlation function may be written 


as 


wT 2 


R(t) © Sk) 1? Wa(k) Jj{4zek sin S-|d°k (168) 


Substituting (5.1), (5.2), and ($.36) in (163) and evaluat- 


ing at t= 0 


Se as 








eee ~§n202k? 
2 x x 


eo i Gide fe Sn dk 
Be ik ener tke) 
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sin? (7k h) Pa vk 


o,28 f= : dk. (5.31) 
b (mk)? go (nk, )? y 


It can be shown (see Appendix F) that an analytic expres- 


sion is 


a2 


BR 4aBo gk (w5a) (38) 


R(w,0) = + 1-Erfe [22 pecs 
a 


ata? 


Cae 2 ou vO 
+ as le Erfcjac + 25 | 
a 


+e % Exfe [ao - z| - 2F rfc (a0) | (5.22) 


and Erfc is the complimentary error function 


Eric(x) 


The computer program used to solve equations (5.6), 
(S.7) and (5.9) and calculate the infinite summations is 


listed as Appendix H. 





ae eee 


bg ' set of parameters is’: 


Sb i 0 meee ee ee EO 


Numerical results constitute solving for the signal and 


VI. NUMERICAL RESULTS 


background coefficients and using these coefficients to 
specify the form of the filter to be used in the threshold 
detector as well as to calculate the probability of detection 
as a function of signal-to-noise ratio (SNR), background-to- 
noise ratio (BNR) and number of nutations (M) used to make 
the decision. The probability of detection is determined 
for the averaged equivalent signal-to-noise ratio and for 


Specific point targets. 


Many different values for system parametcrs and target 


locati 


ons were used ta compute the coefficients. <A typical 


ane 


1 


nutation radius = 15 
width of detector = 30 
height of detector = 1 
blur circle standard deviation = .5 
background correlation length x-direction = 720 
background correlation length v-direction = 720 
position of detector in detector coordinate system 
x direction = 15 
y direction = 0 
target coordinates (if desired) 
x direction = 0 
y direction = 0 
pointing error standard deviation = 7.5 


A detector utilizing the parameters above would trace out 


the area in the image plane shown in Figure 6.1. 


} 


All units are in milliradians. 
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Figure 6.1. Scanned Area of Nutating Detector, 


This type scan was chosen because it crossed the center 
of the coordinate system which would be the center of the 
target blur circle for zero pwinting error. Also, in the 
case of a random pointing error in the absence of bias zero 
mean is the center of the coordinates. This is important 
because the threshold detector design is based on a Gaussian 
pointing error with zero mean. 

A sample calculation of 121 coefficients using the para- 


meters above is listed in Appendix G. The equations shown 


in Part V were used to check che accuracy of the coefficients. 


Table 6.1 lists typical values for the summation of 121 coef- 
ficients for various sets of parameters and the calculated 


summation values. 
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TABLE 6.1 
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Calculated Sum Actual Sum_ 
Signal (point target) -6826894 -6826993 
-6826895 »6826994 
~5438002 »~5438087 
Signal (averaged) -2990664D-6 ~2991232D-6 
~2328138D-4 ~2328142N-4 
.2651595N-1 ~265159Sp-1 
Background .6013477 ~6012230 
~ 2347420 ~2347420 
~8335525 ~8319764 





Besides the gencral shapes of the background and signal 
spectra, variations in the spectra due to detector size, 
background correlation lengths and coordinates of the detec-. 
tor are of interest. Curves in the following figures show 
the effects of these variations. 

Figures 6.2 through 6.5 show the envelope of the one 


side of the two-sided background power spectral density 


a~ 


Sa (w) = z 6,6 (w-kw,) 


where We is normalized to 27. The curves of Figures 6.2 and 
6.3 illustrate the effect of varying the background correla- 
tion lengths. In Figure 6.2 the correlation lengths are the 
same in both x- and y-directions. Practical measurements, 

however, indicate that correlation lengths in the horizontal 
(x-) direction may be longer than in the vertical (y-) direc- 


tion. The curves in Figure 6,3 are for this case. 


abe tae 


ar: wake scart. 





E 
The curves of Figure 6.4 and 6.5 relate the spectrum to E 
changes in size of the detector. In Figure 6.4, only the 3 
detector with BNR equal to 9 and number of nutations cqual iu 
to 100 are shown in Figures 6.11] through 6.13. Each figure 
shows the envelope of one side of a double sided spectrum. : a 


Two cases are considered. First, threshold detection systems 
designed for different Gaussian pointing errors but with the 
same Size detector are shown in Figures 6.11 and 6.12 and 
second, detection systems designed for the same printing 
error but utilizing different size detectors are shown in 
Figures 6.12 and 6.13. 

The curves shown in Figures 6.14 through 6.16 show the 
matched filter to a particular point target located within 
the scan. Notice each of the filters exhibits a band pass 

characteristic. This is because of the additional noise 
component (background). The oscillations observed in the 
point target spectra of Pirdres 6.6 to 6.8 have been sup- 
pressed to provide a look at realizable filters. 

The probability of detection for the threshold detectcr 


was shown in Part IV to be 


] be 


Q, = SNR VM} r |f iss eee 
1 + MQ, BNR” 


xl? 
k 
Figures 6.17 through 6.20 show Q4 as a function of SNR, 
BNR and M. 
The probabilities of detection for point targets located 
at specific points in the scan have also been plotted. The 


probability of detection for a specific point target is shown as.the 


5a 





rere 





width of the detecter is varied. As width is increased, so 
is the area scanned by the detector. Figure 6.5 shows the 
difference in spectra betwecn a rectangular detector and a 
square detector which have approximatcly the same surface 
area. It should be obvious however that the rectangular 
detector scans much more area per nutation than a square 
detector of identical surface area. 

The envelopes of signal coefficients (magnitude only) 
for a point target are plotted in Figures 6.6 through 6.8. 
Phase information depends only on the location of the target 
relative to the initial point for the nutation cycle which 
is along the x axis. This has not been plotted. Changes in 
the magnitude due to changes in detector size are shown. 

Figures 6.9 and 6.10 show the envelopes of coefficients 
for the signalthat is averaged with respect to a Gaussian 
pointing error. The bandwidth of the averaged signal is 
much Jess than a point target as expected. The figures show 
the spectra for two different standard deviations of point- 
ing error. 


Matched filters (magnitude only) for the threshold 
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Figures 6/21 through 6/24 show Q4 as a function of SNR‘, 
BNR and M. 
The parameters of the nutating optical system have been 


varied to show their effect on Q, and Qa: 
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VIT. CONCLUSIONS 


The problem investigated in this thesis was to determine 
an optimum processor to be used in detecting targets with an 
inetared hutating system. The nature of the system divided 
the problem into two areas: optimizing a spatial processor 
according to size and shape and optimizing aw temporal pro- 
cessor that takes the output time varying voltage and decides 
if a target is actually present. 

General equations for the output from a nutating detector 
were known from Samueclsson's work, however these cquations 
had not been solved for a specific detector. At the same 
time it was observed that the fouvm ef the covariance fun 
for the noise offered an easy solution to integral equation 
of the Karhunen-Loecve expansion which made statistical 
detection theory appealing. Some work had been done on the 
form of temporal processor using statistical detection theory 
but this was limited to only one nutation because of the 
common background noise component between nutations. 

Harger's derivations provided a means for describing a 
temporal processor using statistical detection thecry that 
based its decisions on multiple observations and ingluded 
the background noise. His work was extended here to include 
unknown parameters, amplitude and position, which more 


closely characterized the system. This extension, called a 


threshold detector, wis shown to be eptinua an the case of 


~~ em 
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small signal-to-noise ratios where detection is most diffi- P 


cult. The derivation also used tne Gaussian assumption to 


“2 Coptobent’ 


describe the background. This may or may not be correct. 
But in practice, the detector's performance may not suffer 
greatly if this assumption is wrong. 

To implement the threshol’ detector, the integrals 
describing the detecor output were solved using the Gaussian 
quadrature method of numerical integration. Checking the 
computed coefficients by a summation provided a means of 
accepting the validity of the integrations. The form of the 
threshold detection system using a rectangular detector was 
determined and the frequency spectrum of the optimum filter 
was shown. To specify performance, the probability of 
detection is nlotted against signal-to-noise and background- 


to-noise ratios and the number or nutations on wiich the 


decision was based. The signal spectrim for a point target 
and background power spectral densit,, 15 also plotted. 


A designer may find these results useral in developing a 


wae 


compared to the optimum. He may also use these result 


it 


I 

| 

bE System or investiguting the performance: Of an actual system 5s 
] deterinine the size of a rectangular spatial filter to be used. 
Me Future work shculd include extending these results to 


spatial filters that are circular or e:liptic. The sane 


numerical integration algorithm could prebably be used. The 


threshold detector was found te he highly sensitive to back- 


i 


i 
processer shouid be included te mininize the effect or 


ground correlation lengths. Thus, some form of adaptive 


tN 





pe ees 


TT I I gem 
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background. At the present time, a measurement program is 
underway at Naval Weapons Center, China Lake, Califernia 

to determine average backgruund correlation lengths tut for 
different environments, these correlation lengths could be 


expected to vary considerably. 
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APPENDIX A 


DERIVATION OF SIGNAL SPRCTRAL COPFFICIFATS 


The radiance on the detector is periodic because of the 


nutation therefore it can be expressed in a Fourier series, 


a jnuot 
H(t) = EH e QA) 
n 
where 
T -jnw_t 
1 © 
Ret | H(t)e dt (2A) 
But, 
= fxr 2 
H(t) = fx'(r)e(r-e(t)) a?r (3A) 


Substituting (3A) in (2A) 


-jnw_t 


+ 
ey dt z(r-c(t) e : (443 
: ee 


4 
A a aes 
He Ja rN'(r) 5 


The transmittance, 1, may be expressed in terms of the 


inverse Fourier transform, assuming an infinite image plane, 


; -j27k- (r-e(t)) 
Meo) =jaerst@e =o e, (SA) 
Rewriting the exponent, 
k-(r-e(t)) = key - pk cos. ot - ok sinwt (@A) 


Fquation (4A) can be rewritten, 








T -j2 ker 
2. achieiel a ese: 
Hoe dy N(x) rf dt farkrt (hye 
: j2nAsin(ut+9) -jnw ot 
i “e € 
| where 
é 2 
| A V kx + ky 
-l 
@= tan ky/k,. (7A) 
Observe that 
2 : eT. Ys 
nicky) = fa’r N'(e (BA) 
thus 
. T -jn(u t+e) 
< 2 mA SeTeVerR? 1 1 ° 
H, fark x Gare (ae rd dt e 
jnAsin(upir®) 
‘e (9A) 
Using the rclation 
: i eae a f eine ei tsine ds (10A} 
; ne y = ae 
° 
the expression for the coefficients becomes 
| 2 j Ong. 
= ey cack) etn Imp ko eke 1 
3 HL = fa2k Nt CR)TA(K) Jy [80 ye ee (11A) 
Let 
2 = 
¢ = tan kv /kx 
then 
gg dh 3 eins 


Rewriting (PIA) 





: of + 
HA 2 ;" fore N'(k) 1* (k) e LA [ere ky tky |. (12A) 


SO ie atte tn ssctarsbteet abies cate te, ee 
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APPEND] X B 


DERIVATION OF BACKGROUND CORRELATION COEFFICIENTS 


The nth power Spectrum component of the background at 
frequency Nw, is given by 


l T jnw, 
Ry (t)e dt (1B) 


Y 


averaged with respect to the starting time. 


‘ 
Rg(t) = Af efpcesrpacey]at 
} | J 


Using the expression for B(t} from Part Il, 


T 
Re (+) = 7] dt pfs cde (rye dd ay 


; 
fs Cy) atry-2ep] ee (88) 


Interchanging the order of expectation and integration, 


| 
But R,(t) is the correlation function of detector cutput 


Rp(t) = Af dt joey farts a[ry-o (ter) | 


. ‘s Edye ye 
tfryce (ELS es 2] 
Assuming stationary background, 


[Scr 8" Cry) * egCry-ry) 
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where ¢p (4,75) is the background covariance function on the 
image plane. 

It is now convenient to express the integrand in terms 
of two-dimensional spatial Fourier transform, To begin with 


j2nk'r 


j 
t(r) = fr(ke 7 7 a? 


k (6B) 


where t(k) is the Fourier transform of the aperture function 


(1) and k = (ky shy). Then 
t[ry-ecter)]t [ry -e (td) 
7 j2n(k,'r,-k,'1r,) 
= farkyr Gy) fereyrt(zye 72 22 27 


-j2njk, o(t+1)-k50(7) 
catnfiy oC 42007] - 


Substituting (58) and (7B) in (4B) and interchanging the 


order of integration, 


Ry) = fa?ky vy) fark, eA CK) 


(88) 


The last integral becomes, with R = ryor 


: -j2zk,°R “jecksry 
[fetes cae 7, “| =~ (9B) 


But the bracket above is the Fourier transform of pit)» or 





the Wiener spectrum of background on the image plane. Let 
joe Rae 
WC) = fre Tes et (108) 
92 





ieee 


If BOOK) is the Fourier transform of the point spread func- 
tion and Wy Ck) is the Wiener spectrum of background on the 
object planc, then 

ree : fe ry 

WECK) = LEGG TS WOK) (118) 


This relation is analogous to the output noise spectr a 
expression du to a noise input to a linear system. In terms 


of We (kK) the last line of (8B) is 
jen (k,-k,)¢r 
way fo 72 7° “2 a?ry = wey ky) (248) 


Substituting this into (8B) and performing the integration 


with respect to k 


Ra (t) = fa’x [tk [* W§ CK) 


fT darks fe(ter)-9(t)) es 
oF e€ dt (13B) 
0 
Since p(t) = (p cof.ot, psinuoti, 
. -o(t)] = K w (ttt) + xn osir (t+7 
k [o(trt) e(t)] Ke cos as j aye Aw tt ) 
“k,pcosuyt : kyesinwt 
= A(t)cosu t + BCv}sinu t = cl(r)cos(e_t- v(t) 
3 Qo 9 
where 
A(t) = P[k cosu 1 + kys inet : kd 
Roy os -k sinew t + kocos._t - k 
Bit) cf ky sins, hyces 5 ky 
x Uioipa mo RCo Baek led, co ciety payne a oes 
Coij = VO (1\+h (7) = Oo eae = COSw at) ChY KY) 
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and 


W(t) = tan’! vey (14B) 


Therefore the last integral of (13B) is 


-T j2me cos(u tt¥) 
J e dt © Jo (2nc) (15B) 
0 


—j— 


where oD is the zeroth order Bessel function. Therefore, 


the correlation function Rpt) is given by 
ae et eat Os Va 
Ra(t) = firco] Wp (k) Jg(4nok sin “sa k (168) 
where 
wo tance eo, fe ab ab 
k yk tky and d k dk dk, ‘ 
Finu.ly, the power spectrum coefficient BY is, from (1B) 
= Ve weed k 
B= flral? wyaoa"k 
T wit -jne_t 
. 1 ‘ ° , 9 4. 
ih Jg( 20 sin “| e di (17B) 
where a = 27kp. The last integral becomes with x = wot/e 


T 
J J (a Sin x) cos2nx dx 
° 


7 ae ME Ee 


aH 


T 
- ij r) Jo(2a sin x)sin Int dx (18B) 


The cosine integral equals J? (a) and the sine integral 


vanishes [2]. Therefore on using (18B) 


} 
| 
| 
= 








easier ans 


k-plane. 


2 teeme wee tee eee tows cae ree meee eee 


- 2. 44 
B= fire t? [Fock L? We(kyI2 (2mpk)d7k (198) 


This is the general power spectrum expression for «a nutating 
detector output. The power at frequency hwy is expressed in 
terms of aperture transform t(k), optical transfer function 
FiCk), Wiener spectrum of background Ww, (k) and nutation 


2 ‘ ; : : : 
radius through Ji (21pk). The integration is over the entire 
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! APFENDIX C 
EXPECTATION OF hy 


t 
| 
The essential calculation required to find 





4 
Ay, C245 = Patty. (2,433) (1C) 3 
bt 8 
i for Gaussian background B is the expectation 8 
2 
T Te 5 - 
I = Ep exo A(t)B(t)dt - c f Br(t)dt * 
° ° ‘ig 
: 
| 
' where bs 
; 
t c = M/N 1 
and 
$ 
aE | 
= ” 7 
A(t) =] 2 (2, (t) - $)(t)) (2C) 
n=l] 
B(t) may be expanded in its Karhunen-Loeve (K-L) repre- i 
f i 
sentation provided the mean square value ; 
2 3 
E(R°(t)) < @ (3C) 
and thus the integral equation 
do, (t) = ; Ry(tsu) ¢;(u)du Get 7 (4C) 
; can be solved. A; are called the eigenvalues of the ecuation 


and 9; (t) are called the eigenfunctions. 


Using Mercer's theorem, the correiation function of the 


background may be expressed as 





ge eg case oe: LS - 
R, (tu) = Ano, CE) o*) Cu) (5C) 


The K-L expansion of B(t} is 


B(t) © E bye, {t) 


with 
t 
be = Sf Bt) oF (t)dt O<t<T (6C) 
5 : 
where the cocfficients, bps are uncorrelated and because i 


B(t) is Gaussian the coefficients are statistically inde- 
pendent with means We and variances Aye 
Also expanding A(t), 
w 


Alt). 2-5 ae (8) 
k k*k 


with 


TNE cree AP ene gether ot RoE Ke 
ry chi al baal Me a J dee ut 
2 as BS ae aay 


T 


a, = f A(t) ¢#(t)dt, )<t«< T. (7C) 
k 0 k ae 


Rewriting in terms of the expansions 


Te 29 chee aco 
1 é 1 
I=E expr ys (5 °£ ta be + = 2D a*b.) oe. (tyo®(t) de 
B | |: 2a k 2 k Sack 2 k k n 
JT 2 @ 
-¢c¢ f £e bbe ¢ (t)en(t) dt ‘ (8C) 
o nk > 


So acon oe Ro aes Oe eee ee 


Using 





eo we ee 2 een nev + erent TE 
: A CR AR IT ER = ee OE es re ee oC 


. Ty a 
to) J ‘ 
1S tee jes 7 (j abs + 5 akb,) +c : sara | (9C) 
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Using the independence of the coefficients | 











OE Te ot year te eae. | 5 
Co J Eon CXP E a,he + 5 arb cj), 
oo Ja \2 any)? 
= exp Lyn OxP -c |b, - | (10€) 
For Gaussian raisdom variables \ 
exp[w m m&/(1-2wa) 
E(wxx") = Se Wh 
where 
w= F(x) eS a. Se “ 
ho EL(x-my) (x*-m$)] (11C) 


and w is a constant. 


Let x = by : a,/ec 
and Moe oe a,/ec s 
then : 
a.,2 exp {cl -a_/2c]°(it2e\ 
Eun exo | -<|%, * 3¢| | : —s “ta § 
= (1 + 2cd,) a 
Substituting (12C) into (10C), and separating terms 7 
i 
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Pp we ae ee 
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ar > 
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eT TIE er SR 





p -1 ee n u 
J 2 (1 + 2cr0) exp 1} : rat" | 


Han “a n/®) 
C « 
exp ° > ae sy =i 
, n 
wu (uk-at/c) 
mee | “5 7 “Tee: (13C) 


This result can be written in closed form by using the 


definition Ay e (45%) 


2 
1 la, | My N T 1 T d 14€ 
pe 13 2cK ae May Satay to) ee 
where 
eye Ne 8 aie ees (18C) 
Qty a*2 Non 2 + 2caL gels aye "2 
js such that operating on a(ty sts) with 
TT 
Ry (1, ots) (-)dtj dt, (13C) 





LSU ReCeyats): 
00 
and using Mercer's tncorem yielc : 
N : a 
oH qa(t,,t,) + : Rg (T_,t3) q(t,,t,)dt, % Rp (T),t3) (17C) 
Likewise, 


gr = ; + i RoygK Sy 
ee an/c) Ky Os af/e) 


1 + 2¢2 
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tyo 


m(t) - ACW cay at (18C) 
cq 
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O's me 
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Metenttietie | 


cartel vibe tii tildes ae 


= eteee 
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eryee anes +: 


1 RRB! STW. fon 8 
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— toe 7 2 TL Eee neg, een 
where 
- au, 
h(t) = = Tee ¢,(t) (19C) 
n n 
such that 
T 
Nh(t) + 2M i Ry(t,.t) h(t,)dt, = Mm, Ct) (20C) 


Taking the logarithm of both sides of cquation (13C) and 


using the equations above 


eo 
tle (1 + 2ea.)7! 
n h 
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n ,! 2 2 
+ Hy f dty 5 (2, (ty) : s(t 


Xx 


M . 
aa (Z(t 4) ~ sqt2))| a(t, ,t3) 


T M 
N {2 
- f m(t) - 54s © (2,Ct) - S,(t))pph(edde 
- | b ac gel 2 2 } 


int © £ (l¢ 2b.) 
n 


» 71 M M M M 
+ DL be pital ke eta jo Sb. Seti n ck, Fests) 
NH 0. [een * 2 maa ™ 2° eg BP ey OP 
M M M M 
- © §(t5) F Z2(t,)+* © S)(t,) © S_ (t,), q(t,,t,Jdt,dt 
ej ogee ee ge hs ge ae Pe de 
T 
- J mi(t)n(t)at 
° 
7 
pot Fee ty SS, (2) Cet & (21C) 
M R=) . 
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The likelihood ratio desired is 






2 M T 
A, ((2,}) sexpyy - if Z (tS (t)dt 
: n=] o 
1 M T 2 | 
-=- F sf So(tjdth ol. (22C) 
N n®] o B 


Taking the logarithm of both sides and substituting the 





expression for %, (1), = 
2 M 7 
gn Ay ({2,)) os eo 4 2 (t)S,(t) dt 
M T 2 
-=- F&F f S,(t)dt 
n=] 0 
= 2M 1 
—- \ 
+ : Rn (1 + N k) 
} T 1 M M 
+ If | S Le (tay a eas) :# 
NM °0 kel] aed m=) ™ 2 7 
Z 
M M M M + 
e Yo Se(ts). 2 2 (ts) So oSatyy cf) beh 
; ge] Rg.) m=) ™ 2 a 2 gel & 


M 
S_(t)) t Sn fts) q(ty.t,)dej dt, 
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£ 
£ m=] 
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ge > a M Ja 
- m (t)h(t)dt + q © (2, {(t)- S,(t})h(tjdt 
° b mI Oo k=l : ‘ 


(23C) 


: | 
| i 

| 
10} - 





eres tee See OS SS eee Se PE 


APPENDIX D 


KARHUNEN-LOLEVE REPRESENTATION OF EQUIVALENT 
STGNAL- 10-NOTSE RATIO . 


er © ee 


The variance of the threshold statistic 


M T : 
We ¢ J Zz (t)q(tjdt 
n=1 0 
where 
: ke Ts 
a(t) = = ie - J f(u)q(u,t)du (1D) 
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under hypotheses Ho and Il) is 
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NM f? -2 ee eR 
var» MY @Papar os fa-act) faud wk, (tw) (2D) 


Substituting (1D) into the first term of the variance 


expression 
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Using the Karhunen-Loeve representation 
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i 
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bye TL ty (tog (eat, 
then ‘ 2M 2M X 2 
NM -2 Pee ee a Wk eck 
k ck 
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~My pte : (SD) 
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The second term may be expressed in a similar manner to 





give 
a ait Boe 1 ok os Grae: 1 i 
ne a Ad (Hare (8) b fei{ are 
N “k Nk 
fee eoeNo we Se 1 
VAR * Sp Lf, ff ——— . (6D) 
k 1+ Fr Ay 


The equivalent signal to noise ratio in Part II is 


To 

a : f(t)q(t)dt 

des ey (7D) 
vVAR 


The K-L representation of the numerator is 


M = i - . M 5 oy. if | 4. aN ) 
roa f(rv)q(t}dt = Laser [f, ‘(1+ = ss) . (8D 
n=l n fe) ayy k k N k 
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d may now be represented by dividing (8D) by the square root 
of (6D) 
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APPENDIX FE: 4 

; 

NON-ZERO TLRMS FOR NUMERTCAL INTEGRATION ; 
Uquetion ($.$) for the radiance function for a point : 


target muy be written us 


-2 2.032 
a dof," ry 
] + J 


ae 2: sees 
*¢ jn¢ Jnle \!X +y" ) avey 


er sinax sinBy 
x y 


where 
a § w/¥2o 
ft = h; Y~2a 
c ® ¥ip/o 
fo 
ps i (te Xo) 
v> 
Qe o (Ty %o! 
and . 


a (QE) 


¢ = tan ~ y/x 
Let E(x,y) be that part of the integral tnat is even about 


both planar axis which is 


n Dk es ; re 
aa hoe “X°-Y° Sinax sinky Poke \ a 
E(x,y) a ray J, (ey +y ). GR 


Using Euler's relation the integral becomes 





I= [E(x,y) feos (rx#Qy) +jsin(pPx+Qy) llcosnd - jsin n¢)dxdy (3) 


Using trigonometric identities 





rT; dxdy 


i=l] 2 


b= fecy) 
where ' 


Ty = cos Px cos Qy cos n¢ 


T, = -sinPx sin Qy cos n¢ 

T, = -sin Px cos Qy sin n¢ 

Ty = -cos Px sin Qy sin ng 

Ye = -j Sin Px cos Qy cos né 

Te = -}3 cos Px sinh Qy cos nf : 

Ty = -j cos Px cos Qy sin n¢ : 

TS = +j) sin Px sin Qy sin ng (4E) 
{ 


To save computation time it 135 necessary to determine 
which of the ciaght integrals are always cero and for which 1 


they ave zero. 


Observe that 


4 n 
SR gt eR ie?) 


cos n¢ 





Likewise 
sin no = (= he I (x+jy)” ° (GE) 
x“+y m 


The first term in brackets is always even and may be made 
part of E(x,y). The second term in brackets may be expanded 


by the binomial expansion with the general term 


2 pe ee te 
ee ee ee ee s en tee oe Ce ae 
e be oad i‘ A 


Gr? Gir s (78) 
The first integral of (4E) becomes 


1, = feCy)cos Px cos Qy R, (xtjy)” dxdy (8E) 


——- iS 


For Ty to be nen-zcro, Re (x+jy)” must be even in x and even 


Will always be zero. 


in y. From (7E) this is true only of even n; for odd n, I | on 
Examining the other seven integrals vields | 
t 


fECGy) [cos Px cos Qy cos nz 


+ j sin Px sin Qy sin nojdxdy n even a 


JEGGy) Leos Px sin Qy sin no 
+ j sin Px cos Qy cos n¢])dxdy n odd 
(9E) 
Cbserving that symmetry about the x- and y-axis is necessary 
for this integral to be non-zero and integrating only when 
the integral is non-zero permits one to integrate only over ia 


the first quadrant. 
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The signal for a point target was shown in Part II to be 





H(wgt) = [N'(r)t(r-p(t)) d’r . (IF) : 


Assuming an infinite image plane, H(w t) may be written in 


terms of the Fouricr transform as 


eaem OY 39 


H(w,t) = Jxraort(Ke a? (2F) 
For a rectangular detector (2F) may be written as 
: * 2 2 pe aes 
i H(w.t) ae ese Bye sintwkx sinchky : 
QO T x T y 
i 
: > ; ; » 2 “ t 
een eer pda Re hyo! d°k (3F) 
Let : 
x = 72 omkx 
y = ¥2 otky 


and noting the integral may be separated into twoparts 


ooo a . 
as: x(pcosu t stg X4)) 


2 
1 “Xx -. WX 
H(w t) = = fe sin —— dx 
= a J V20 
sin = JI : 
2 2 j= y i te “y 
1 oY ee 29 y(esinw, (ry ure 
TT y : 
(4F) 
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Since the intcgration is -~ to +” only even integrands will 


be non-zero. Rewriting (4F) 











WH : 
= 2 sin 2 is a 
2 2 -Xx G i 7 : : = 
H(w,t) = rl e = cos (2 x[pcosu t (r, X9)1) ax 4 
Fe 2 sin a is i * 
5 e BS eases 0: ve , ieee 
= J e cos . yfosinw t (ry Y9)}) ay. . 
(SF) 
The general form of the integral is _ 
R = R(B) = a a she cos(Cu) du. (6F) - 
Differentiating with respect to B , 
anc) , 27 4-v" a 
di = Si e cos By cos Cu du. (7F) | ¥ 
oO : 
But the solution to (7F) is known [2]: | % 
ae e/a le te 7 (BF: | 2 
Since the constant of integration R(0) = 0, 4 
a 1 dvee)" di gee a 
r(B) = +-|f e ev+f oe dw] (9E) 
F 2 T 0 Oo ’ co 
Let a 
v-c iy 
Tae | § 
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and 











then 
B-C Bee 
2 2 ¢ 2 
R(B) = 4 2 f oP gp + 4 f eo 2 ag 
“} mw J-C/2 vn J-C/2 
ar Ui BSC Ase g Bee 
where 
2 
2 A ey \ 
Erf{(z) = eS dy (11¥) 
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Substituting the result (10F) into (SF) 


H(w t) = i Erf Lee epcosw t - r_+x + 1 ; 
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¥ 20 
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The mean square voltage of the background may also be 


solved in a like manner. 


From Part V, the mean s juare voltage was shown to be 


ef = firdl? Leoa!? Weck) ak (13F) 


which for a rectangular detector becomes 
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Letting 
t = 2nk 


the gencral form of the integral to be solved is 


= ede 
R(u) = 1 i sin? (t/2) eo \ 
me 2G t2/4 ae +t 


Taking the first and second derivatives one has 


BZ 
dRQu) 1 {- sin(2ut) e°% e at 
du T J, t at! 


and 


m 
a? Riu) . 2f" ¢ cos cos(2ut) ont 
du 70 ra +t° 


The solution to the second derivative is known [2]: 
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2 o7a* ; 
oRt) = 5 ie Erfe (ac- 4-) 
a3 a 25 


* e% Erfe (ac + 4) | 


where 
Erfc(z) = zfs 
v1 Zz 


is the complimentary error function. 


hoo 
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t 


mentee talib sok: Dalchat ibaa a ae, 


(14F) 


(15F) 


(16F) 


(17F) 


(18F) 


(19K) 


With Jengthy and tedious calculations it can be shown 


by using the relation 
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v [ Vv 
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with initial conditions equal to zcro and 
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The mean square voltage thus becones } . 
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| 
j -_ 
112 i 5 


SSE i 


TY a er 


— Smee eee ss 


APPENDIX G - CALCULATED COEFFICIENTS 


DOVVOOOOT MOGOVWHOVIAOGOODNHNO ODVAYSDVIOVEIAINGVOOVOO 


~ 
Bg er 6 8. LOO ee eee Gob eh ee eRe O. ROe CATS ew OL A, Bela rele ig ire 


‘= DODVTIOOVOOCHVVHOG sOACHAVOETOOSNVOCOHVAOROOO 


i 


m 


_ 


LOVE MED COVES SP TIVO DON A|DOI MN MOMS ADO OKORQHN 
SO DOOD DO OO et a et et ed et et et AOC OY 
(rrbrveer ta dte rears reseaetaererrraeroegedaed 
badd WIOIITIOOOC HOM COACHINOOOMQGUHOOAACOOMOOOAG 
OO DOH OW QMCR OD FOO MWOM NAG FOO SNM DONMOMMOMm NEN 
CNet Dt OID WOO SD it CO ST tet FO IND EEO DOR OOM WS he 
DO DAT UN DD OP KO SMOOVE ARR OAD NAD OS ETO OLE 
DO SF eAaAGlOPOAMNM VOT OM NGI SF Dat ND OG OM OLS DQ me OU 
FUOAOO ONE MOO Got a OOM OINOG™ NOANNOODQ AD NOS UNV et 
PAD AAMAS OD APF OMAAIAIIVGMMON OO DAL D— DORMANT ROMO 
DOTA PFOPN FIOM IGN DNR MNAAONDMAAS SNOSARONANEHOS SM 
AINA OY met ER OOD EEN 07 ed eet BG et UN SF RAM FOR NIM MRL AIBA D CAA EE 
ere eter ee tee ete eteaev Xt o> Fe *eoeee eee teeeee_oee te 
PEED OA SIN SDN QO eo O EMO AN OMIM Cat ED mee ORNICOUN 
\ ‘ () ' ' ' 1 ' ‘ { t ( ' ‘ ' ‘ ' ' t i] 


oe 


~ 
- 


tro | VBPCOSOOS0H 
wo @eeesvwe 
ext [MOOD IVUIG 


me 


OOWDVONDAGINOVAONVODO0OAO 
Moe he a a ee ea ae 


DAAIANMAIDIOCBOVOWIOOOS 


10 
Q 
‘). Q 
g 


) 


IPOQDDBDOAAIONVCAISVIOIALOAIAOAES SC ONSTOBDBVIOOOVAS 
phavorvevrrr yr rerrarr terre evar ta ea 
GOOCH OOIC O25 BIO COUIOO OOO SCONKESCOMNCCCe 
SQW VENER OVAL R OIG yh caer DIA CNA Se mmitio re 
“oo ‘Sea ee PPHRTCKAAATAee soos Da DOWD 
LMAO D™ Dam DTaIHLPI DANNY OMS OQCAMAPDAM a 
, ror SAS Aw ORR MAD YE ~ ap gine can 4 mF VCS AN. se 
LS ea NEM RTE MAO N DAD N SRA DP MONAT NOE 
Daw DUNO AK QAM ARTS OLD OL ARNT AP CRAM Aa 
PLAN FP ONO DQ OOS PNOINN D PaO INN ee Oh OMOC LS 
OOD NOCD DODGER RO DO VOPR MOEN USE Oa MOm Eta 
esoeoee*spteoeevt* eee epeexeeeenere geeet*eet® env eoeeee*ese 
ttt PR PALO DOM REM ODO WMNMA SAO ete 1 
( { ' i) ' 4 ‘ ‘ ' ' ' i} ' ' ‘ ' ' ' ' 


NI BIN UN NGA STOTT AL SEN ES OT FR AEROS AR a OO 
1 


2-92 


ie 


eo 
2 
° 
th 
~y 
a 
v 


om 


T.0914,0980-92 
GL 5% 


-1.06 


AMUN MMAMAMMM So Seta PUNK nnnnnwn 
ODODODODDOO OVO OVCSCSCAOGDSIDIaSSAdOCHQ0COO 3SCDAD 
trb’_beetpprype tary pyER GH dt eeaeda 
Om oaewcardcelenragsaec SEOG SS CeCOSeSask & 
HADDAD OW MONA O AN AADA AR a AOANAWID LAG 
OOM AON T OAONOMD Dern SNUVE WO AMIE DOIN 2p 
Nap of NER DONNA SE OD eet Nee er seem 
COte Wo OT HON me. OP Pi peta am Nee oar st TMARS DS 
Belen oO OR RMR OMNAC EON Se MDE MOON NS VA SQARe Te SOM oO 
n MOR AAD MOAN ONIN GHA DMNA TADARE RA DOS AOS Fackmar 
am OT IA FOU MM CONAT RR EON QANIAANOOONNNGD OLN 
OeWAD CIDA SHA Py AN OtDNNOT Piya DOWER Fenove 
eoeeeneeeoeeteee ese tee ee eee ett © te Ces *# ee Fee ee © 8 6 o 


AN ON SO He IO De OT A et ett I OLAS SM QQ at et PS OIA 


n-04 
59-04 


(eo 325-08 
oMeEL-95 


- 


EDA SOR ZSHO-4CL NT NOR MRO ANM PAO OPAMYUALNOE SD 
Ed th ed tT OCU IAQ INO CR ONIN SO BET Tt <A Nee 





OQDOAVWDDOINDIAONNOTHOOCAGVOADONAL YyIOODVVHBIONODONO DV dDQ0000F 
S$etce@se tee Fe eee tee XH vee Fe eeete etn teat ee Gaetan esas genase 


QOGOMQOOOCE DAIHDOBHOAAT OC COSTHVOOVDOOUNDNHAVBVOAMODOOGN0O0 


PS TIN OM DADO aN AUN FLIONIN OOM OO PO tC TULA OM OP OO MA ERR O 
NN INN AINA EE SOT AA MMMM MAM MMe MOMS ETS et re UAMniNnguh : 
Prvruve dg evrbrbteeterarverteran Prbtrtarpteterepturreedreyre iad 

ODOC OA GOL ORR AEN OOOO AMO OUIOAMOMOQAMOMIOOCCIDO a 
DRONA FP OATS SAMO TOA OE OO OIVN DS POOR PSNNO MOH gy Pw 
QMADAND PADD QO FEMA DDG WO DOF MNS ANP aA MOOIHAPA SE ROD GN YEO 
DWFDN DINAMO WAI a DD PDE PATS ee PIN O SF Det Qo PEN DNS HO DOR 4 
wd PRT ONE MMM DO FP ROS CHOIR OM QOH OQOMMAMN ONS POM ROMO 
QSL FE OD POT SRY FOUND OMORORNUMORANRES EU SPOR AN COGOM 
TAN SND Sete Meret ae GODO SE ER BINGO RON DWONST OOM DARA OOM DO FUR) 
DAO OA SAAR DA PF ODF et ON INDOOR MN APR OA Le HAANIMNOOAN FROM 
DES EO RIA CORR AO EIN IAI OIA PID ced at 0 me SAID mt et FOR (Pie ny 
@e tee tee eet ee Feo terete even tee tee Fe teeteo et eee egeeeeege 

DAS OR AP eS EO MRI AE NO MENT BIO DO HMMOOMNA EPONA HOY 
be #€ Fo thrrvtretrrrprrreepgererretrrverebrebtrretoretya ' 


oe eG 
a 





cmoet ie ues. lsh 





is he 


iolelslolelelolslelelslololelalelelelwielslelelololelollelolwlolslotelelololo lel ~lelaimlatnlalo) 


eeea tee te tee tee Fee ee Fee Feo evat eet este epee voe ne Hy Oeee 


SMVOIVIIOONID AHI OO VDOO OHO AD VOOG ON OVOS DON AVDABDGADANDO DOCGO 


er cla ts 


SPADES AMR MA MOM ATID DM OAD & Sor 
BOON ONG  -ODWDOONV OS IDAVDIGOOAQ 
: ety bt pete ter tebrr eet torr sverae)d ‘ tharyyppretorres 
BOG LOCO MOAMOQOMOOLOAAIMAMIG PSOANGCE FOOLS AOL OS 
OP HOS FAYE NN TNE DUN OM MAY EC TR Na atm ot et ett rR WIM 
POW AIS TRIO OOH KO aD SOT OS Fe OO YIM NOD FOR CAN the, 
DOI NAO OrIN ayn eA atts DAR ON GUM SSO OD tm ot eet tet Ot 20907297 
DOR DNS PN EO te A Det I CIM DONA AE RPO CA ARNE ERAS TOT 
WO HO DMI OOS TIS MORO SEN FOAM ANRANR G O ND SS CO RAIN Ro men 9 
NOM WIM OMNS AED OP EE DORON SORE SRD SED RN RO MAU ee : 
MDM PENNY TNE NASI OM US EM Va ORAM TT OM SD OMA DONATO 7 
DODNAIRD AA FOAVIA MOY ST OAM OMAN Ram og OF OAM Rey : 
ef*eeeteeveea tee eevee tee Ge teet ee teeteeat*e eee eenvesee © eg 6 a 
ENN VE POM DD LON SN NN et Dt EN OS mt NN NN et et et tt EI mete IAIN oe 
t :) | i { ' t to ' ' i] ‘ ' ' ‘ ' ‘ ‘ ‘ ' 


THOMPSON PLAWNADACI 5 
SGDOMIANIOVOOIA79OD9 


-04 


wou 
we 
OL 


0-04 
~04 


rn 
9-04 


rf 


es Ob ke red me Be 
og et a! de a 


C& 
8 
oh thi met pet lade kei eee he pm et 


REAP EE OCOD OOS 
MOVOVAVNVSOOC5SO00O 


~ RIE ROQVWQOV00 ON attest ee 
o 

pitrprreptdetespet oede 

x 

>) 

vw 


Ooo 
¢ DOOD DDG ag eh rh ed at red et tg tt 
4 pearee eee seeepgerrbeetebarva 
COC OODDNODOOCAEADSTCOOCU NL. 
De D SF DAL D OCH OLDAMSAODMAMNDIMNMSNOROEA Hae CHA AAIN 
BRANNON ST ON AD SOW IAO MAMI DDOM £ OT NROIND RON G OD KAN 
ROOKRN ADV CIM DING FCNNVE SA OTFONMAON OND NGIONIUIM F ELOMANO NGM 
POD Ng mt Be SAID ADA. O TOMO SOY OMNI D OM AMANO TIRE ODIO PMH Ne 
MIE NOM MOA CL SWI OOD NAY FOC IE MTOR OF AO OM LOD MPA RE O-WCaAN 
MOTE DTODAR MIG DANG CIDADE Oe DRAM EASE OME IYO D DMAMTNOHOR GO 
WO ADNAN ADOM MS TE DHAVIMOM DIO OMAN FH DARE NO Sete3S DO IO 
POA FRONON ONNMVOSN DROW NRONBSROFPROFOOTOrE me ORR te f 
eevee ete eee ee tee to evs et ee eeereerrere ee esvnvrVxuevesneerseeteereeeee 


PON NA tt HO GD TINA SOM OS OY ON DOIN AINA RM TN et EO EIR 


0& 
%-99 


won 

oo 

Uk ' tt 
ODOIOQOOOAANECAOMOCCOOCOdIOGEe 
OP NWM AE = 





WD DED DOME ADEs DAO MINI NOR OT ISAM EOE DDO NIM SAC : 
SPT TOMOANAARU SO OL DOOD OG ORE REAR SAE AR ROW OD aww D 





Oana 
wevys 


redithwe. eal ot 
tee oo lee ee 


114 


; 
i 






















DOANCIDOVUESAC: F®ODDODONOHNAGEGMOOOCSOVO0ON 
e © @ © @ Be Fe © e wp Cee Fee Fee e emt Fete e_eeeoee 


QOADVCAOOC SBVVWONO OOMDADOAGDEHAQMOVCCOMCSO 


ONS 

=n 

tod 

O20 

weed UL) ood 

orm 

ome 

Naw 

QUetin 

Qor 

Nt ey 

MIEN FORAIOADADAASOAAAODAAAATDAIOOADACANQGOAM 
efeee_eoeeeet*®eeeteet*Qoeut®Fenvneteaoevpt#eetQeevr ee @ 
SON eee ese oer or eee Core ere ese . 
' 


cet ee pe ee 


@ey @@ @& @0 © Fe wee F @ 


ier wee “ee 


DOO 9DBVBO9OOCOMDOGO9900N0900900 
° e 


WOO LODE EOD DI~ OUR ODYMUD YR Or CRE RR 
MNODQNANIYVOADEI AN DO OOD O DO VE AVBDONDGOBANDNND 
: ; yt a A AA ak DA Oe Ar 
; BONOG ACU OMADOOCIAOOOAEQOIAMOCGNO2G 
PASO aly ey HNO FOOT TOLCKROMIR OONN ONO 


: Oat AO OTM $F OT LOO AEDMA SHEVA MO 
‘ DAW DO MPYOD ONE M OOD ANTON 4s NOR Ar ae 
i DOT FOU SFO SINOMMAMBAN F HOR TPH LPOOWPAND 
‘ NOOK SH QO CON OMAOMO WNC AMAR FARO NR We 
WADAMO 2 IDNR ND OMAN Fame BS DOO 
UN COR BOD NAW AES DE AVN CDNAS EPMA DLA TRAN TRY 
ODA GN MOND SPO SR APOANE SOAS AM CAD MANDON 
; @eevwet®e#eeeet ee Fe *eeetet @eeer*®ekteete-eete 
' NARADA MAD AIA Ma MANOS A IAT ON Mt Pe 
t 1 ' yt ( perde t ' ' ' t ' 














+ tetera ener ese eealem LE RERD s2tee 1pm te onde be heme on 



























ANN CAM AGM EG MUNA MOS PSPS PEPE TUNA 
TTT LEE LELOERALE e oaee Gan oe 

ti ' rhorrerre yer oer ed yea rung 
AadGAAAdAdAddd>AdAGEdAdASADEAL AAADAG 
DOBMOMN ST FONSMMMAAKNONS™ DOD AOATMORD 
OAD NU PFANVE OVUM US OP ANOFDOMNINDODROM 
MOALFWIDOIANINGMIASMUINDOFNRTORROGMD aes 
ROMO IONDNOTVOMNN AU ea T MONA AMT AUIS NA 
MOMGAWTVCE O FNAMAM DOS THPOGMEUO DEN OW 
WND OM DOING DE LOR DO Ne CDONO OF OMIM O OAIN IN 
FOO DORND AE ODAORANENSROMRONDEEOODDM 
POO MEP RMU TON OD SIN NININ CPO SIMON OOIN E O 
eaeetarvteteevermtteevreeoetrtetoteoeeepwet ere seeeoe 


ON tt OED NIN i dt ed md tm DO OUD POTN mt ett PLN TO 






AOAN nrnokn PFOaNA TADS IT OKANO OF OPO 
DG poe SPP: FDQV OG COO DO San tt at ad eet ad 


mg me me md ed ed rt eg OE ed net Oe gy ot od 


© 
@ 


nh eet attach ehti Line de balla, cnthe ast mmr om 


PON aa tL eRe pega ne im 


APPENDIY® 


3 


WY ~ 
ere ax 
Zu ws ~— a 
wor On Zz uw 
wnt ot wa 99 Wh 
wiw co om Ou 
hm - ee ul 
uw Of OW Wy oO 
ere Zo aus cod VU 
wig D uom =“ an a 
Onmr.g Nw it an es ~ Zane 
QIK Z ~~ oO wo ud Ds 
Qa “ND -—e of [oo Le) 
a Oe Le tmwU pe ad 4 w Wat 
waa wood ys — ww) Z2Z2aocn0 <a 
se ww Onw CO> w v NOR WwwRe 
wi lla Camera amr ae LX 
Dus.t wWweUO Yo 7Zui tr amt 
Coney Ze yk eas22 ZrOOeZNALEZ 
tie b=? TUEN>SOW Z2ZTSIOS DUS ee 
CMs ee er OF ODD LL tone LQOO>O 
ut De dm wae ee Ste et aa 
xzs><-7 CORMING Skee ee OO Ke 
milo le 2IZr~ Do SOU Zid Pa ZZ 
D wT NM SWwuldepev ern BID ee 
weds a 7h CMa lm mei 
ret oe DU — 4 oCsaiuwn Ios 
a-re Zoemowl 49024 Pirerurme DD 
arn Tew thee OD D1 OH ee addtr bOU 
Se ms tie <> 7 Zeldeuoy 
WS ah tec wae ew Te et Zihin- dd 
ob Ss as oS ee ROCK ICV 
ate wie CS meet slaare st 
Ott Wh UFR CDZ 2427790" aa 
Cou ius (G2ZD Aint ite CO Ser aH 
Der (COW Be eS TU IN ere 
mou 2 eis > ws wig VEaEE 
dem Cl Come Tete YS LOTR eR 
De sine Pov we avon IOoitrh 2s 
“s ORSZowe D NOSSO C rd Ube 
fe Zl aeTI Uys CU oe wat 
CA Gee Der eT LT ak ‘sofvda CoC 22 
cea—- Clo tense Se es CCrra yr se 
wee Uo: t lei Le O SLA Lye 
QAwLwt< Ce lA mw IZ Tees LOCC AEN 
ine ¢ “2EOS27arenw dete Vb Ff OYss 
wld Coe S ASR CUS Sib bheeK LES 
wb -"o QZ Lers Paw ZARLLoOoOw lt!’ 
eke o> a od cH WM Mw bem VdCGaace —C 
ATdJsc° » 3 we ZF was 
ae oe wt wets u pe A ee oe 
=~ Orr +O >Sor>ce 4 =>> 
CnZow aZ22714aT=- ~ucoox swe da 
eCduc dduutuevislagarxtX exec ZZ 






ae 


a 


ee a ae ne eee LO EON Ree. 


oO 
» 
- QD 
o a 
~< ° 
o 2 
» ~ 
Ow - 
De Qo 
>< ~< 
“iv eo 
woe ~~ 
eul oc 
—is e 
Oe < 
sw « 
mi ° 
—e w 
2t m 
> 7 ° 
Tan bas -~ = 
“~~ ux<e o oO a “2 
t oe we e Pd ° m~L 
cm Oc> ~~ ’ roo c 
- Fed < 9D K © Gk? ‘5-1 
Cc ke \ 24 4 —-“— OFC Fr 
lmeee ST 95 - =m mz er. Qaim 
ZOrwOr CC el2aZ Dm 240° O09 ud ~ 
@NOT IY Ft >t ~ ays we emu oO 
SANDOR TET mm eu ke OZ © ° 
owe Ce ee Qe ABA eN Ya DUO 0 
Ads EN PDE aeoS ame 
whe IIe eH OI we Hu ame LOD OCS 
Meet teeMTITNan awl ate wrcim sm 4 


beg Q409 or Zhe xrz ~uier oe we NS eD 


HE Se >rZe se vewwtrt2Za agwa me ot 
SroanNAAIOVON SR sic Frm LU Zero 
Say RNAS OK SUR TED oe eS SPRCALANG 


ore ne © 28 POUV II Yew bee www wm 
CYMZFZLCNAAM LOC UNNEM avo>uin Re EE 
m2 wT OCS She wHead wutdaaa 


wero 


wh SAF TOC Oe wR wer Ea 

Q¥E RSV F ade Zh Zoya ewwe ee KH La OAC 
Z—-CCecOCutidx<reaqdidarauucarytCOltocreZd 
mMOOVOCOLeE aa Les XO Been GCF BU ee vw 


ANS 


> 
~- 
ew 





SUBROUTINE BKORNIVP,V) 


GS 

Z 

fe a 

> oO 
ee OE - 
Duia’ YU 
rat a 
wo > 
uO 
ca> z 
on x 
wiz 2 
FOO om ~- 
Zn -_ iu 
(Baw > > a 
- Cc, 2. ww 
wa 7, 2D S53 
me Oo qs Zz 
a ww > <a 
uae - 2 
wDD =z 2eF 
oO a W]~ ft 
ozzy or a 
Giy wea 
Lue) we ws 
=e Lae 
Revo =m OW 
Za COO0 Ff &- 
wD ° < nm 2 
Wi eu rou wm w 
~- oe CG oS _ 
52.7 “wo asc YW 
aGMnu DUN 
Lease uD u 
Or ws an eek Ok 
ud - f2zoO0U0rtw 
—>~ ax0q raga 
Fude area ZOKU 
Carn Oo Ouse 
San Lusso > of 
xOWNZ Oran Cw 
Lisa FZ9ZD 

Ube NomMuTIas 
woomM Wax FeO 
2rzx2e Devedew wt 
—Douw So 
Reo = sr x 
Du re >C ult Spiny 
OOuDd = IN 
ax 6 See 
Gow @ ><NJ s 
zS<—-= & d-qad 
Mate NN ZS Om > 


Ae NONE Hoge pe: 


pHeALe bE eRXeR Vy X0q YO 
c 
U 


) 
(50) 
25) 


h 


2) 
20 
2 
H¢ 
oP 
wl 


Lm—pe ws eo mc e 
CeO e>OCEN eCO 
PON ie OMe OTS 
DTA~ se nUaw~Z © *O 
t~O. WN eo oe 2OO © 
alfn~FOoCcunaro 
WO e ZDIT OZan oN 

em eZXA N HOn~ 
oOnmOn Nawal e 
2 AMNSNNNEZ OO 
Sew EET FY e0O 0D 
Soe mri ir SOOO 
iw>OoAraeaGesrOmM & 
x e eO 


COS (-1.00) 


ee ZZEUGIO NM 
Shae ESSE Serr i 
QtzZ20UTEz=Rr oom 
Ferme OVOOL OCU 
—COGCVVVVECOCOszZa 
mn 


QOOVOS VYVYU OU OUI SOU 


CALCULATE CONSTANTS FCR INTEGRAL 


- 
a 
~ 
_ 
a 
~ 
N 
s 
” 
ce 
_ 
wy ~AN 
‘apd Sd 
wo oe 


t QO 3WwO 
71 esi OT 
HSN iwwre 


PARAMETERS 


FOR SPECIFIC 


SET ANNULT 


VWVYIVU 


Ord 


210 


Nm 
af 


v 
mc) 


G 


a 
if 
mx 
mil 
NUS, 
—_ © 
ae 
=O 
Ne 
2 ~ 
“Nt el 
muiNRt MN 
om me 
NINN ret 7 
moe | N+e>D 
Zaconnm2 


TeMAXR GI TC 210 


‘Su! 


iow WNoOo7~—~" 


“tm 


maw = 


Ow 


Sah N EN 
aww 


LNFOwKG 


ON OZON—O 


150 


oO 
Co] 
“ 


~_ 

-_ 

= 

Zz 

=z 

< 

ua 

Oo 

a = 

ures 

oOo «a 

=z ™ 

> 

oc Ff. 
+ 

i oe 
” 

x x 

a <1 

=z = 
2 
~_ 
N 


VUUVY VO 


wn 
2 
~ 
> 
< 
=z 
a 
= 
v 
NI a 
oa w 
iu a. 
= Oo 
Oo aw 
- 
Vv ” 
S ue 
co - 
” a) 
< Zz 
- <a 
ue 
a a 
- Zz 
ay <1 
da 
a rea) 
oe - 
a = 
~- ws 
"T= — 
ot Lu 
(oo ¢ 
_— Oe 
Sym oes 
om > am 
cC~— eo Oo 
New a 
an Y 
~ ome 
res uw 
Oo aad 
--—O 2 
zany aq 
mt = 
acm Cy 
aaa Y 
VVVY 


117 





a: ae : 


a a eae 
SRNR ER 9 ene eter m. 





(MB hae 


sey 


~ 
> 
2 x 
be 
z 
vy bf 
a VU 
“p> ~ aw _~ 
~ w” ~ br -~ x m 
ao > uJ wz ~ a VU 4 
a - ~ ~ _ = Zz e 
z oO x eo ik & - z OQ 
< > z - ~ ” vo N 
=z a pr x _ v= "n ~ =“ oO 
- - - zu 2 Zz — - = CoO 
= T Zz zz ” _—— ~ ao a a * a 
pe oO 5 CGC ~~ az > “ z= «@ ~ o 
- <a oo < x >i Zz - = - oO ~e 
fod uu OK uw - > zZ x -9o x< ui an 
- a ~ Wy <a <= a> a te ~ Aa 
-_ a >4 aw Ww wes = z x Zz oO oe 
~ oO _ a QO oa au =— ~ - —_— _— — - ° oo- 
+ aw uy ate ioe Y wh vo « 7 QE Mes -~ WN NA 
=“ _- ~ oO ~e a < “ -7F Zao a4 Load msm IOEN 
~ zx ue ra ~ FZ uw = _ Zz> C7} oO ~ Use on 
nw zo oO wrens et can) — DZ+ -_ ~ = Ce Uns. ow a NOK «K © 
- 20 -_— Owns De - COx 2 = x - >—-s ~ + aae = > mae 800 
-~ Zo ere Ayxr zm Ke ow a O aim 2 om << -—O —>> = D> « OnmWMNO o& NS 
_~ OQ -O O Zwei Zz oO PS tee ee aes 7 > Owe = Be Aee => zo 5 ~~ eon NSTC 
~ om = «§ SyYa2 wn o— IN aw enw D> Orr Oe Ze ZAG “A © O DD NOOK Re 
md ony 2D mayen > aD Zao— | e Tor SoH ~-Q we ZZ - oe SF mm mtt De ow Ht 
~ oz wb 4*MNT Oe wiuw 2Z"o wiuwetOeZ ewa eoNteoT lee ua WN 8& OyrmUxXKKS 
< ~ Ze ~Nawwro 7+ — uw =a Vow MePODe I> ew >ONsrw & SDE PEOS SDNADOWS 
< a rd uJ u,a NO us Z7uas ~ Me TM St Dh AE VDA AAS a ae ewww eww 
=z - - WNZKO™ « - Ovoeoez 2 <a Oc G O en Ceonnwl— a al O4D2 Zeierkerrere 
Vimat io DPD © MUNZO > recom erm > Of>7UR-0rrF COLOKGZ-“ObFr— D SPOZNrk Hts 
J om mt QO Noro 1 OO QU eae t eS S220 TOA ZCOrZ2ZF VU aM DrreSDZOILZSAaYH aE 
a ——_— ~~ = woh ~Y = VO wwZ = Wd >Re Oe Swe Lew twet 2 4 SF ZESe eRe Oeereae rudd 
a LKOZ Q OdA>acau oO OwVOZeO > ©& rloOocwMTOCnO#waAuOoragcgxy Lf DBO DVAX™WOGCOVUIea0!Z% 
oO zo00 YO CumMwor>« Y OY >Ccooun MY OP7> >A ZOAQOS FPOT>RV>raAGdyH UO NMORRMVGCuUbLHUUKEUiW 
oo oS 
a oO oO ono oO o oO oO 
Col “w ~ MASS wv oS oO oO Naso 
~ ~ ~ - 0 wy t mI NOT tt 


VW OVvvOV VVoOYV vOvovVuU WY 





te 
Pere ene cee 
UG ree part a neem a 
4 
' 


t 


<I 
c> 
= oJ 
—_ 
za 
re 
woo 
wa 
mo 
re 
MoD 
ur 
wd ~ 
za > 
owt * 
i —D 
DO <~— 
Us moon 
ar —/o 
2 ~> 
: > “Nia 
: xy om To. 
‘ ft Us ~ eer ny 
‘ ~—m wn eZee 
, ae ' a-7 
~ co — > 
Lt er Loe 
=O De wha 
<a be Noww 
- <0 OAM 
7a ~~ + om 
~ <ac IT nan 
>» a qQ aQ's 
wn /~ v. 
x= oi 25 ae 
—_— < ar DO¢ 
D> Bs ue COW 
uw ubk e vu o Om 
U xa CoO™ 
Zz ZY -—~ - ety 
oO Ourm mm Ourllcr 
- —t— to VWUZOFT2Z4a?7 
- = Die mls 0 © wammtiy 
w wua=s= JIF mK ea D 
< Zea atewwstr O 
=> Doe MOta uw Dus 
wm 2euUwW OG rem Co 


VOVOVU 





AT eee, cero rm 
RH ng 
tRNA Zen eg ye MTN oem ree me cee oe one 


{ 
| 





CAREY eer EET ON 





re 


Te 2 ee eee Ne Clem aR. 
; bg an Le MERE EE AEN ARON GF PREM A OTT 








aw ° 
oO ° 
, — Ne 
' “n Oo +o 
b a N oO 
; z —_ oOo 
7 iw a ee 
' -_ N Wo 
W<et - «ad 
_ -~ Cc “e+ 
f wa °o > co 
, Lew w - on 
i> ~ oO oo 
5 oo _ ~N < eo 
wo ah . ~ a2 
wn > a & ~ —N 
} Of zum Qo xa - +O =. 
wz ar SS N ee eo folola) ™ 
“iw a 4 ~ uw Cc zo Sa 
XL D2 ny - © « oo t 
2uU Fe 2e Oo wWM e «%O at) wy 
oe 2. ee w > w woo a c 
wu w at w ~ e« oO Oy as + th! 
a mw - =z © - Qee <a N - 
Uru) wae wo Ovo Coo ut ~ w oO 
xro Ue he » W NOD oO . Fr ~ 
eu Ler wy ~m—-o oe 200 uw wo < oO ~ 
- mam Lu Ovte ure en e oO _ 
oO mcry COG Ue NNR od HOODS Zz “ «a N NO > 
+ WwW Lid <a Ow — aA~-OAKHO OND. _ _ a Cre i) 
‘i uw ro eu Ure N eIN «= © oO enN™ - a - es 
' of Daw Go Oo ve ££ Oi «xX TO re a a Ww r- —O 
: - Cua we auu O Wem >oar e2ON oOo og xrad _ cS 
% w Fw CMNOuWw wom aT en OY? NOOO. m— Ww Na Cart Fs a oO +-_— 
4 w ora i > oO XL DT nor W802. O Os. +41 _ a « 
+ a udu a OO § ZN>e we YOO te O wn Wig Kx«KOO re) ~~ ~xX 
i —_ ACOUO qt Ow oe FOG. COBYoOn -_ - 20 oagnNN ww us Ww at 
: a sr “Od roe wm Oe FAD AZIZ) eG Ht OD z Zw bie a  lx«K NIL oN 
; < LOZ “er ZOw Om es Zar at vat «Om to < an xxOO “i Oo <a me MN 
I Zz zen Vv cwo aNO—~ N&O my - oo - aro “~O = wo 
“2 Ore LW Ee Do RANLTONSNNEN 8OO7 Ort uw CG a sow a Ne oe xe yy 
' _ eZ OFDM _O Sew NUTT SVN PO ee IO z. —~e ~ aanNN oO aN FW Fee + 
‘ ww itm zee SZeen7ederwsCOO's owns oO aN _ Qno~-0 a v- ONO OY Om 
e vw Vv Ome UIE WONG y GQLAmwGOOM si> AN ® Ww ow a. e oF-e LL nN Ome wi te 
; ul wit w2ZxFDexe w O>rReur Ze eOnNn SOO or ~ Cay ta “_ —-— mANR ON ON 
z <r Toe oe Z2 2a GRENIODOL KUN ws we aw Ww. C7 Oy BS Oe ay eae . 
_ -— 6 =e Zz Own ROD, SNS NANO XK OIA OM mm WD Om «Oo © euiad Dd Geawrw NOS meas uie 
c red a0 w pte se aAoOorqnodaa a rnd dao WIC 2 © AQ etl ee lee 
> 2.2Zz- >CHuH ot CNMI ZTLZSPSN tCZE- EO O 4A W>Gviwnit OO>>e4 *¢ Z GON 2 00m 2” 
Oo Can = i ReZZ2ZYIUITOUUANO NE MO Ft D wet NVR eK edd KK < Wit Om ANDO Zee 
a =a” to aa oO SUWGLWISS ES EFREaAnNonaone x OQ Been ernn 1 bende wae NINE NE ZS 
oa oa i XxNGCU QVe sees SS Vekr mor ms ON 3a £— FTZNntiwwtt i tt - mt ae ee 
~ > <a Q aww Sm MeOCUOCOGUVQOQOGU Yi at ALONNO KKUUKX u“ NO~ UNX Owe Oe 
“” May NN=E Soa =OOOVCOVOUUGOCOLZ2aa VW Oe NEDOVG Ori BOW WY OON RRO ONR UG 
! mN 
i oO ° 
\ w oO 
i ~ N 
i VYUGVVUVOUVVIVVOVY YOvVoYV VUVOe 
iii eres Ua ech ult 
Pa peg ee oo nae a NR i a ll 





eo x 1 RR rpc py earn Men eine apie aon: 


i Tier wei Dm ke re RD STE OR, Cat OI ne ame ot Poet 
Fatt y=! 
. | 
} : 
F | 
a | 
Pa F 
1 
= 
4 
a “ 
E 2 
i 5 
4 5 
z. Zz 
j < z 
Ey go a. 4 
of . - 
4 ! 2 
‘ ! 7 < 
: | oO w 
> 
j ow 5 
° i - 
x & nas 
\ td o 
x uw @ 2 ios 
c uw S s 2 NOS 
cs zie > < oO 
- a = 2 a = = Te 
s . (are so > + + 
© oO xz cz = = = SSS 
- Ze WwW ND a4 = 3 = 
Ss <i oe o > a x g Vit wero Se oe 
~ as « ° vce ua - co < ee ee 
_ ox Wh o -2 sa 7 o 3 = 
= ao ee woe - o w zz a f+ a—-ny ¢ z= 
Zz ok. moo Fe Ge © ae an x HO oS 
> eae oO + Oo us w eZ + = VU . e “ee 
<— iw — ~ oO _- ~ Tro ~~ 2 < Pee "OD 
on 1 we = OQ -~~ Oo ze SOre- Zoe us WwOO> mS 
ue (te zN > OC: Cy = O ~~ swesr 2 oe oe OOF SUN) 
Oo OK - - “OOoC¢9o =~ wr Tt — Zz Azn<—Z eo > x mt et mA Pre 
pesee fou? nas “D* a x a e mR +r Aemvy Zt mw oO aQabeOOF7FZ rier ew 
co mt om -_ ee On O72 ~oceH~ o WO > weNaeww wove aw ST we 
we | Ge) B+ AZO wenane = 2 ¥ arZau 775 u code ER Oe 
oa iNew ans eon nH tote OMe D> RB meme Hm aoe Oar OFZ NG 
= ONIN Geh kf S& OF Cty HOS GO wzrZurrZies i BOON 4G 8 000 ~ & PAHO © 
ee ee <a 7 Bem L Oe Re O KZOnnZa cfu DBD BCOOm ~ear20  eOT2Z sGOL9 
S oN $e wo anugn Zee YM eA S WOK) ae weet ID at eee KNW Tet OO 8D 0 
+ O & BRA Nh eZYOLNK eZ uU F&F LONORZZUZU CH ZZ 1) OOOOdM~~Oa Se a 
tt -CmrEOo > Wr LNW TDI DWM EZ Zu eutrrsetr > OtateutOrsroi Oo Oo 
MROZTO Of Chee YN NCO SG fF Oe few ce eb e mH OOOCCINVNZELO wen 
<— <x Ke Sere ON rl tO SZ Owe MUU w th DES SS = eww Ome 2 7 5 
SZ BOGnxOG OF SuaerOfnevsauget UO Bae OZUWOZEAZGOO DD NANO Gbunword ~GCITN 
= FO22AG OO UZZOOSG-Saswemoaa OC CAOdGMUnKAGASGSLS YF >>>O2MRraazZaaaG>>>>> 
oOo =aNWO w ACs S 
mt A yds ? 
NS ANNAN co ee * 
VU LOGYU VEeV"Y"N Sr 





++ ena eee ee | 





a 
Zz 
~ 
> 
a az 
z 7a 
~ eu 
~_ mS 
Aww 
c7nZza> 
eZ eZ 
YAO 
—+wem 
L723 
Orang 4 
+>4>> % 
~sunNus > -~ 
Farm + a 
wernono 5 z 
APN OA o oO 
> > aA > « ao 
~~ ewe zo oO 
mO-Or MW = ~ Qa 
COOOYU »SZaw~ Oo 46 0 
Gmee +002 QO We ° 
© 0Om aUZO _ > oO 
O *O em Di wo KErO nV 
Ohm eh eee MeO Fu 
x Kuwos xKCcwd oa Da 
a js Meree sOe> O«7z—~ Ne 
= om om” aw fe “wD a 
o momo TVA OoOZ0* —o 
— aZzaz2 ow7Z WHC HZ 
Hew eZ eO Gu ee mQix bd oe: 
mH eM ewe EOS CE 
Ae Ie Dew eet et OO OHO 
O04 ~ Iw 90700. 2z2Oa >w * 


OX Z ZZ PF AMK7 AR RR mUOY AY Ue 
OA GODODOR~e ZZ RrODwn 1 ZS 

Be wwrwre TUM ZU ON ae 
O§Huvnwd Guwg cCOOCmMOwu DL 
OMS we JRO COUP OAH & 


Oo ow °o 
oO ON o 
0 we @ 


REC(T) 
OFPQ 
1) SUM=0.00 


UM, SUMJ 


6 


I 
§ 


NS,NUM 
+ SUM,SUMJ 


) 
, 
2 
’ 


WIPROWOZIAZ+ SEUWO 
ates Fe IO 
w 9900 VUDUNZ 
OWNER CRO LNW ONZ eww mero anda dade daa 
PLZ ZDROZ NTP S Zn SSE ZOD TT EF LLS SS ways 
PeHSaKePFywswL VCS lLitZeare CeCe rCeCCarvousrra 
OOHKDAWONA3LDOUWU CA DADS CW CDISDOIIWSOQRODOOZ 
OBAMA LCOPAVMm MOP AAVOMVA G COU tee bbe LUN eeULULW 
_ 


I 
I 
9 
+ 
8 


( 
1 
1 


CALCULATE SUM FOR INTEGRATION CHECK 


950 
285 


15.6) 


Ulew 
NO 


—ae 
oOnw 


— @¢ om 


Ge2XeT ae 2Xel PEELS 262% 


—~™ @ 


atu rs 
WO OQUIRUID ees 
MH AOUIN ee eo 
ODOUR AO KK KE Ke WK OE 
D mt DON AON SAUD ILA CG DOIN 
CO ee ee ree ee nee a ew Ne nee Oe 


RR eR RR ee Re Nek 


NOS OOO 
ANOS AO tate ett 


*,°15.3) 


122 








RBar tenes camer en emeeee ine. 


FUNCTION F (Ry THe NN) 


BYVVVOLU 


~ 

wr 

nu 

— 

Oo ” 

ao “os 

a é -~ 

- ~ 

-_ “e: 

+ —D 

~N <— 

~~ wn 

oO I~ ~ 

<= —_—> oO 

—t s~ +) -_—— eC “—— 

ox —o -~—— ° =“ 
ee -a <~~ 77 Yo > 
ml re we —— o —— 
' aps< <a =r a xT 
oO —a ~~) Ree <a ae 
- Usw Zo woe —_ aon 
x Ov mmo oO ZZ x aia 
' ae w Vw a 22 4 wa 
a Oo aa aw — ow a — 
_ me oO weetOoQg O NZ oe od 

rd te amo e—- One Oo mes 
Oo mn -— aQas er UN Oo Viti oo 
a eG O ”aReenn mHreuns OC OC DWH ao eG 
ao wu OO MK>—K>e |HBMNOD VU + ODO. @s,rnr @ 
Suraan owmw 2-85 O9OC7C « anmm Qmom © 
ue ae hoz moDTrKKx™Ms*«K SO? KM HIND te 7990 < 
xo mn QZRracndG «¢ «KKM Comw AO ow 

Mame ZUINS ewww ww OO OO eet we 

we. COON 2 ODOR ANZ LN. e@sTO!} PWe naw KSG WK 
me ROO eee SOOWGK Ht th i DIDI M DH %Matl Dl 
ULUNZZACO MUU NMNOVU DY Zeman Z 6470 ZF 44207 
Kec te eZ CODOCO oo * NIDWK MOAMTAZAOWeA es 
SBaAWwFr ad —-OCOZ Ae 7 ut YA xX eww we MINA ORUZDIAIOFUD 


QEETTR Mee LEK KKyH Ie wwwool~- 


vn 2Ow- 
TOR Ot TOUS HNMVVRMK ULC BDO GAMO VEWOaNo nw! 
mIOOVO PG OwOxX ODODE Rea OA DUOC DUE MOG DONCw 


o 
~ 


wuZCroa 


oO oO 
Lama) + 


see eee 


2k eee ee He 


SUBROUTINE WAACR1eR2¢Re THe WeN) 





a 
Zz 
<< ye 
wo 
nw 2 
ee Oo 
xa 
Ow zz 
p= <= 
Ww a 
zD> - 
© 
“na uu 
-> x 
aw Qa 
~“ = 
Ous 
a0 a 
— 
al =z 
In" “~ 
> oO 
os a. 
<I 
mo a) 
Zz oa ’ 
wd ww iS 
=e an 
n-O wD oo 
w za 
Me Zw —_ 
wy ACcCoOoO fe 
me — usu 
> MOM o 
aus Mie DIDS 
ake = oO 
rR LMU) 
a2 Or ow 
Px 2 ” 
qu mt be 
to OoCrr2Zz 
z= a ow 
z ww MU > 
w Wise Ww 
2H xrawcrlg 
- “Leocze 
am" Od aqwre 
Cr) axe Fw 
Lm a uo @& 
ano * re) 
<2 a xX uw 
Vig 2zarstZe 


COCEUCEEECLU LEY 


-™ 
-_ 
Az 
a 
1 
~ > 
o (oa 
_— oO — 
=——4 e = 
vw mt ~ 
a ee a eo o-_ 
Or wm QoQ —~ 
o omy ca a 
Laem mm uw 
Penn o 2 
due ste = + 
~ £0 mn 90 a — 
oee~n DO EB vo ~ x 
wee OAGCe — ae oh 
Jato) 7 es O Cc. we, 
dete e—-nu vac ean -— 
WaeNes | OF~2ZW _- Ze acm 
a wih SwO rt © te BAN A NON 
ZUENY «eA I emt e& aru 
FOTN SEARMO HS bwaer «> 9 CO Bw 


ot ACNuewwY FT + = et ND > 
GQNOZ ZL es SG fees St Clee 


ab Faq vwZtina—ae tt ah eons 

Qararteruv2weay we Zt ” nme ZrO 
FRmLO qwnes aACIUO ® OKRK YH KHw~w~ Sous 
=@=OXUOAAMVCrSCZrRrLOeCK Cae SCawY 


o 
c 


ce 
a DG ~ 


Fe ee ae ee EI 6 =. 


Sa eat oe 








mb 


SURROQUTINE GLO2Z4 (XeA+C oD) 


XD 8BHAMISHT BNC” 
ot TOOT NAOM MN AD 
aS DAM NR HA Sto 
mt FIDO WMDMOD PUY 
Nwex<X POOVDOPDIOAMNA 


-~ ade AQOTHONDONN NwaN 
A Im MOD IONMYNOAWL, 
ro we VQAN QOD ino 
z= Taw TFORPRMNE OO 
on aoaeq CQoMmnaerranars 


wD YR REM OMIOP MOY 
PQOVO PO HMODNOIWN 
aX POD Pam VUE CMD 
NSD NOMNDDONNG? Im 
mam NAL AL AOTMNOMO 
sa DOONTOR MAM 
ww TRO SE RASIN OOM 
KT AWA CIM A OM RAI 
Kiem NOT AMI ONnuN 
© 6D ANT AEP ORAS RA 
_——~— oot e®e@eee se 
todd eoerteerer ere Fe rX®ee 


Nw~t DOOQOQDVDODNOCOOR 


24 PC 
' 
x 
’ 


. 
4 

2 
7 
A 


ANXNGROSRODANAN GY 
Pw ALPO ND Mm O.IN TID 
Net ema PO GNNSENAD 
"2 <A PAT UII NPM ONM 


Cine OM AES NOCHADYEMOMI LS UY 
CK mw CPOMMM MA LH OTTO aa 
FPN AAR IN ASD DDI co 
CO~ SINR AACIONN OSM D 
Sex -ADQDD GASTRON SAG ++ 
DO Km IP DON COM EOD 

oh KITES PAAR MOY em 
ZL mw MAMN AN TIAP NC DOM ot met ome emt 
Dit X eNO™MOIIAL MH FORK MOO —~_ 


meme KOM TRAM SICHIOI + x Mw 
OA ommaeNSNOODVDONLOD mM AI 
em mat PROC E KEATS EAS via. 
CU Ow RPOTAONAAIDOLO4 ® NRO ¢U 
Wi mee BOWING Ne PE OOO EEO 
Oy mag UE ONE OM OND es COTS 


Ca Wt AMA OKM RK OI NILA NMs | 
Kem NODINMAI MAMAN CD 8 OND i) a 
WILE Meet DNR MOUS NET HO =m nu 2 
a Oe CN FORCE Lown -~ =2 


GAUSSITAN-LEGENDER INTEGRATION ROUTINE 


coc (esses Ci EOD DS eM aa ee 


ONG thai we tw tikes a 


VEQUL 


125 





“oO Tl Sere 


De abe ea AT BE be 
| MRCP MRO e MECN: ot nome wens o : : : PAT Rog NEES Eh * : a a ak ll aI ical all 








Si acl ee 
4 ° 
3 + + 
w ww N 
Z2 2 ~ 
wl -— nn ~ 
Lu DZ _—a ou oO 
x r-O -_—he OF a 
we ue ~2 ‘ 
ww ase nO = Pa 
wa z we wa z 
oO uws< ema . yy a 
oa ad - ~ Vv 
” < cw ™“ mam "2 
Zz Ous wie = ~ 9 Ve 
gg K uwt ~ wo ne 
-e wo a > ~ oO Oo ws 
- Yn ae Ve uk 4 ~ # 9 
wv Zz OY mT w ~~ D cw 
2 oOo —O tw one aQ ir 3 Oo 
> -— eH ~~ ~ awn 
re re Qa xa 2d Oo wn ° ua 
Oo Ze pe em eee hw 2 Oo oO - 
a 22 23> a2Z22Z22u—r-Y - ~ 4 < 
~ iw as2uA ee Se et TL ao oO - Mw 
= WwW mito ww Pofranuwijw ~e ee OD Lo wu 
> vw ee my aed et Zoe ee ee ee a ? 
Zz ue UO aR SURE PF unt a taD DI OO ud 
- wo AUCH A WZZZZ —~™ w& oy (tae FF~— ow oo 
™Ne DIAuY - wa mee ~ —w t1.+ 319 
x ” wWATDHO — G Zz emer IN w NI wee om oe ° 
cd wm Side ee Na <a mem DROW om Zz WUE ARDI Nm Oot - 
oa = me) re (— _ PmOmN Get voy o qd aswirOr*it rt > 
c 2 u -uwd Cu vi Dew 0 0 @ Abe we 8 oc a AOvasiuas wo > 
. a MILI 10 wo 2 ome Re Um ww < ODO) UO a9 
i] = “con Te —t mm IL) ore ~N ws HOO 2. HOT = 
< ia uw MIO tm Q sate eo omen NAGS *-QN0OO y+ ~ 
~ uo CO Site ot ao w ° oe me ee! —N oOo u © © eO>O0r Ma <a 
~ Lepeeeren yD: | ~O ae = -_ Sete ot 2 OOOO Dd 
4 was wer tO oO amt<d - MN MOTAOARAOLON QO 
VT Ws CxS cad “ »~ oO pwO™n ONN ON 'TD wo 
a+ Wire iw A Le om z — mm AISOOADD TI me 
a 6 7D re Poore 24 xc oSaustada. <4 Paes GQ DDD eOKOOLTINA ww 
ac LS 2? Khu Ty CS ous Pig eee elaa aA ™ BR GTZ ein ow 6 ew ow > 
x ee 2ZOWS Wee a 8 (fmm & eONTOSNMN Dm zo > Cm wet TN att OO a 
tw wat za x2 a SO Zany Caaevcoee Oe Vw ef YW Kp Rrww www Ot nw 
2 az SloScec Zz c* WMeweetee Soe eo ob @ wx re sf Pate. “AINA POW x 
- —S Cm-x— rad ee WO SU Lu. dine tu ZOw DOU £2 KOK em 2 
- ela Minuw Su me KK ESC E KY Ot OOOO ua StI SOCwnwtunnNue > cS. 
2 a aIinu CM93 ee ww hee 8 6 8 8 6 oo Om se 1 ec am” Ke Zur ww OrOO!*OD7> a 
3 c= 4 >< KH2O SO Ler rere er ee Eee OS —— > MITC ou *  2e—2 (x 
oa cy We woah 2 ~wLeOdgtice P22Z22722 2/7 2-222 “Se QA NN we Ns wr 
c re ze QW Uwe wer ww ww ww wr ww tt ttt Zz = xo ats ot uw MZN aw 
= _ Canim Zenr LeU KHUN be Ke OUD CSC Ou meade NCOMKC HH “ 
wn WO acacoz/2 COZ mm KZ ee mtr ZTE I OMU YW OnaAddOoO>: vooachacud = 
~ N 
BCeEUCUCSCGELOCEL VoVv YVUUVUY weuuU 


Q-DFLGATUNFL ISX -GFS(NG1L “BF SCNt2) 


r 
eu 


{? 


OHoOan 


SHON 


77 


Tl 
x 
76 
TI 


“ON 
uM apeeree) 
UM BF3C1) 

CN 


COMPUTE NORMALIZING CONSTANT 


SNORM &= Lee 
oo 77 I 
Oe TY ERE RM OFZ (1) 


OCOD ORDER BESSEL FUNCTIONS ARE ODD 


DEAR OE EO.1) BEC TI SBF(TIel-1. DO Al 
) 

a7 _= (3. 0048631) /xX-BF 342) ) &SNORM 
ee 


Beloit ake 





ee ee 


i 
1: 
} 4 
44 





a pal an nen Se ORE ew Tepe 1 St TITIES EE Ae = EMO TCR RET SE BT TE nae oo 


ee Ee Tales Se wre 
‘ a Tee ANE eel 


eb eee 2 enn Le a Art ot ST te. ane 
2 me wee oe aL i 


| a 


1 





roe ee 
oo 290 
AQ O00 
o90 afSd0 
A a 6 oO es e@ eo fo us 
" oco DO So Oo ~ 
“” *tr One oO “Oa > 
us Ort st oD ON ~ 
@a Kece5e fer ODM 83O a 
+r9090 OO OF BO ~ 
a - =: 200 ON NO o 
o ovo OQWnt OM @w-- oO 
we Sao UME aN SO ~ 
MmNITD TOA Me Ae 2 
us WOND LON tre MWe ~— 
z aIrnn AD ND ¢D -~ == - 
~e ew OO Out om OO ~7~ ~|> o 
- ader mee 0-D OW ped -~ NAO ” 
> Ne © err OF AN _ J me oe Ww 
oO 3°09 95D AY ON —Oo -~ -_ Nowe D> 
ax feGO 8240 Ot ND -—O No aw maw _ 
©0900 GOD MY ao <u = Na ~~ <q 
>» mGeet Ore OW nD ~O et quo > 
a WOoVWo COD We Mme Qa 1a wD asQOVY 
a ™ on % Oe FTF ePNW ON OO ae oe ct ‘ard ua 
: -~ _ ~OFTF err WOO 2900 YO ~2at wo: acOo+ = 
: > ou) me eeN Fe ON (OD CO OW “he Aa O+0z a 
‘ Qa Zz CLFOOO 20900 -OFOOF 4 © ~ aQ 4 +ow? <I 
_ «KS I PDADDNATABDDONO DW Om rO at ANON = = 
aq >= = OLHOINDDINCLTOANIA As oO a+ o- OO «x 
* Ie NROANM MAN SOON FOR We - Oo7- +> ACO~ cc m= 
> vu SSO su oO eens Oo + +> Geto us 
' _ 22 CWOSEARDOFSRKRANNER CO BO aS a) _— > awIruU CG 
| < =z5 eK NSP OM DDD 2 oDeoe | 8 Oo + oO §ew Ka PTrVdy= ou 
{ - we TANMMWD OWNS © eN ee —mH O No Ort—QA ame GUUDAVO uw 
4 x " moe e oN ee Oe AOO 09D ON Km —_ Cw om (1) We #,Q > 
; ~ —> AGN eee ee ee QD °GO ASO oO - -~ WW wed ocd rONcCW =< 
' uw « ww 6 #9050 0900 -AMANYD MWA OC --N oO NO44 qa asad e uw 
i J us —~ QDCOCCHOHZONEN © ¢ OD a NS - Ahem ae amy 
= ad WODCODOVGERE BDGOIRNGC a mtty om =z tf - SJ wo. ee —AMNS ‘” 
<i a nC e Fe D ee COD Hehe ww ROR Se mM CO ~o —wCan aN amo» > 
- an AIwSIAS «DBO ePNDOm CY OF a —~ NWN CC wwe wate! co (S'erS ( 
mu Te ee Geir cCNon Ose Co CLF fo fF Ss av - ~~ o 19 (ci~ T 
tu ~ waotrywst erry oP ee Sat YS zZz2ea = ¥Ee > ~ Qa gqco 90c~W~a> - 
z ud ex seer rpoershfeeNeoe OO e332 OWN 2 © Oo Zn oo hemmOQ @ Cr CO AKO ° 
_ za 27TO07F7.C99NT0O09D9D ww BRO Cw HW Cr . OM ~~ Ia 6 me Om oom in ~ 
ee — ZS RCN Socata REI Oe AR ee OK ete DB Clue m-au 1 9 cote Dew PD Pew DOKL ON 
> ro etme Te TNR DNS Se KK ORR Gm te ULL OK ee Ot COD Oe tN OO mw 
G —W CYSRHTAe hs atoOayr Wu SO a thao fe Om § + ed fuwu 6 » ne ee 
w cS mam PV PSN SOON SO ee mm SERCO R Le eC HC>™ ) eet Cw UCU ” veet or 
a 2r ALSO ade SNeOTaIstS pS ee tte he eee EO OE HOU Tr oC 
= aa AER MNS me K— CERT CR Te mw MN Iw He oem >a Ub ypwew > O Sa PADR>AQNkKRK—rry? 
“ =r THAVO NN QAHMOGAOY ANS w WML EN EN TO CON EL te Ct C0 meer iy ou 
ire mC: © 0 0K 0 © OD 8 OO © Oe Ss ee oy Fe Fie et tt CY ee ee I Qe Omer I CIA Cgqeadd{deoaue 


SINO ANO AON AN 


~ aad ec vc 
~_ ~ ~ —e 
COSELE 





SUBROUTINE SUMCK 





— 4 oe ee — ee 2 ate IN ENT OEEAPURUT TI 


2 
a 
_ 
“ x 
A onMm 
uw - 
fod -~ 2 oO ~— 
a ~~ w ~ >> 
=< uo e —— >> /O 
as wm wo oO Cr rwo e 
> = « ~ NOOO OW 
a wu - nA -~aNKNN A 
< - ~ —- Freer OO 
vu ~~ a n" 2 +4 rt - 
— MOY o ~— us ——+ ¢t - 
RrY FEU «x " a2 — «<x " 
Ke WZ a - us ODO eet e~m~ = 
ajar Sua . z ~ — thew D 
Qk wee us w W CC i a a) 
ana aoD ao —_ _ uu JI i—— — 
QW KPe=m ° WO wm Ww JI WM 
~ way <d ~ a « QO 272909 i 
eh pe iD <a we —— Ww + YU “kd pz - vy 
Iu Zwy - uw >> oOo - wx-( tt « ! 
Ru em x Uy oom wo Ft we weru OO 
Cc Buds eo Oo OO-—— + Ub WWw— c 
“nO ai mo ~—a < Oo NN>> Oo 9 OUFA  -¢ ot 
bed ~ NZ e -~ = = Zz _ mY (rer we Mm - 
mu BOW Pew N < QO ££0C—~ >_~_ << ONAN -« 
=O C#S5 cu> acco w Inne on - 2. sQO0 we 
a wk eoeet bee ond Oo —-~—— | a & Cle we _ { : 
ZTE EOT IM cue <t corre Oo + -e TOI Oo 
o> 2<>5 to ow anvoc we ri+in~ ME 2 xeeOrk ie i 
OY nM acre +2aCr W oateawe oOo ae&@Dd ~~ +400 © ' 
~ ~xrI>D2Z2 NZ ¢ > ++09— A an MM Oo9O+%+ Fb : 
Ww hk ue rrd *2AINDD zt aA0>>29 mM tu a «~KOD FF i 
we RreRH @ ~ ~~ <x—~ I cn ~—er 7 
tm ue one OS we KL wwe eda wi<e wh wete~ 
D2 2itz a=eN Onn ae O - -O9N~e9O CO wt O OO4n << 
Vien Lee Agree ONCTO Cocnrc: ees ss Croco « 
aw wel wOoIA9 Om M6 = On20C-YvWN za -M% = ooca e- 
we 2ZeZ Le>Z Vr COO 2 ct *Aree > ox 2 fete ie 
Za mt cad Niko ew WW AANNS Ww “% Ge WA ANNAN XK 
a Cc rN “OO www OI eo ww —~wwow ON 
mw wa - Ouvwow WI wm new wo w Fo > UL www Pw 
Dr eC VC272 nO> +> hme uUbunn ST - 4 22 Fr Uuuu + 
Cre Fw KANC KWAN S> ence rxriw~ DB Feu 3 XCxrxaodt 
x =¢> SWEETEN MIS teow Q Wkuworveg a S2sar QQ whe s 
cy vu QatzLSZEw~w>S~ = Celso = Mer!’ =z OCsSO «a 
TCO Ox ZzOCC #sCdunmu CO wl & te many SC SXkLaec Oo vwuuecuos 
Sie MID =mMRCCSCZECHO- O A4gCCKVaLr we cacy UO §d@cvocuw 
Cc 


~ 
VCC CULCOCOUO CECE CONES WUOU 


Te ee Wee 


FUNCTION R(tie2651G) 


FUNCTION & TS AN AUXILIARY ROUTINE FOR BKSUM 


CCOCY 


{A~H.O-2) 


(-) -60) 


KE AL *S 


vw, 
=O 
=O 


CLOSM 7 Kae KLHMONDNIDNOLS WS 
aI oOw~UuNenwOn 7% 
epawevwlustean 
GQUOUHKN HmdtTH UBS 
Sm MANA AMARANTE NS 
MOM XX >> Mra 


= 


*$16G) 
SIG) 


% 
mOn-O-= 
en INOM 
milk ox 0X 
an anrw~ Nw Nm 
aS mew wm 
mel ath Ayagnd 
Kaew N1OwWDWoW 
wwe: UwCO+C1O 


oMO#STG) #42) 


uw 
t 
oO 
Oo 
e 
~ 
— 
* 
=> 
* 


(PI) 
rY4) 


T 
0 


Os 


Neew ~-O HIPC IOAN 
SaoNATAKNROMO HH OHONN 


Ce tt 
aN aN Se 
feeb Ht OC 


~~ y 





130 


= en 


| eT. 


| REET mere 
- Re eee SO 


10. 


ll. 


ae 


13. 


14. 


BIBLIOGRAPHY 


Abramowitz, M. and Stegun, I., Handbook of Mathematical 
Functions, U.S. Gov. Printing OfTice, T96d, 


Gradsheyn, 1.5. and Ryzhih, I. W., Table of Intcprals 
Series and Vroducts, Academic Press, Y965. 


Harger, R. O., "Signal Sequence Detection Given Noisy, 
Conmon Background Image Sets," TEEF Trans. Acrosp. 
Electron. Sys., Vol. AES-8, pp. 174-185, March 1972. 


HNelstrom, C. W., Statistical Theory of Signal Detection, 
2d ed., Pergamon Press, 1968, 


Hudson, R. D., Jr., Infrared Systems Engineering, pp. 27, ss 


Wiley, 1969. 


Hughes Aircraft Company Report 2341.4/18, Application of 
Detection theory to TR Scanning Systems, by R. Yoshi- 
Tani, October 1970, 


Papoulis, A., Peobability, Random Variables and Stochastic 


Processes, pp. » McGraw-Tfill, 1965. 


Papoulis, A., Systems and Transforms with Applications in = 
Optics, McGraw-Hill, T9080 7 ee 





Peirce, W. He, “Numericai Integration Over the Planar 
Annulus.,”* J, Soe. Indust. Appi. Math. Vols 5,. Nox 
Die GOs? oy June 195: ; 


Saab Report PAD-0-70.33, Description of Sianal and Bac - 
round Noise in IR Svstems with Transtorm ‘letnods 
oes es eae ee pS Nae See Ee dae 
in Swedish), by H. Samuelsson, I[2°C, 





oa Lie ahaa 


Samuclsson, H., “Infrared Systems: I. Expressions for 
Signal and Background [nduced Neise with Space Filters," 
IEEE Trans. Acrosp. Electron. Sys., Vol. AES-?7, 

pp. 27-33, January 19°71. 


Stroud, A. H., Approximate Calculation of Multiple 
Integrals, Prentice-lallT, Ince, T97. 


Stroud, A, and Secrest, D., Gaussian Quadrature Formulas, 
Prentice-tail, Ine., 1965, 


University of Michigan, Willow Run Labs., Rept. 2589-$87-T, 
ASTIA Doc. 460352, “Spactal Prequeney Tiltering,"  ¥ 
by J. ob. Alward, Dee. 1965. 





1S. Van Trees, H. L., Detection, Estimation and Modulation 
Theory, Part l, Wiley, 1566, 


erm ee ee, 





eT a ce a een eR 


ae 
ba 
i 
Al 
e 4 
. 
t 
1 
: 
i 
P 
Hy 
: a 
ae 
“ x . 
aa 
: a) 
a 
. 


132 





